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CHAPTER  I 
INTRODUCTION 

The  radar  cross  section  (RCS)  of  jet  intakes  has  been  extens¬ 
ively  studied.  There  are  several  reasons  that  it  is  of  interest. 

One  is  simply  that  it  is  a  major  element  in  the  radar  cross  section 
of  aircraft  and  must  be  accurately  evaluated  to  estimate  total  air¬ 
craft  RCS.  Another  reason  for  study  is  that  potentially  the  RCS 
of  the  aircraft  could  possibly  be  reduced  if  the  scattering  mech¬ 
anisms  are  well  understood.  A  third  (though  certainly  not  final) 
reason  is  that  many  aircraft  identification  or  classification  tech¬ 
niques  propose  to  use  modulation  of  the  radar  return  imposed  by  the 
aircraft  engine  as  a  significant  identifiable  feature.  It  would 
be  a  questionable  approach  to  establish  such  a  system  on  an  effect 
which  is  not  well  understood,  particularly  in  terms  of  establishing 
the  system's  susceptibility  to  intentional  confusion  or  camouflage. 

This  study  is  somewhat  limited  to  the  region  for  which  Dj'X; 
the  asymptotic  forms  developed  for  higher  frequencies  may  fail  in 
this  region,  and  the  very  low  frequency  techniques  (for  which  little 
or  no  energy  penetrates  the  intake)  are  invalid. 

This  study  develops  two  of  the  practical  problems  involved 
in  calculating  the  RCS.  The  coupling  coefficients  at  the  mouth  of 
the  intake  are  known  in  princple,  but  difficult  to  compute  in  prac¬ 
tice.  The  effect  of  the  engine  structure  is  known  only  generally, 
and  has  been  very  loosely  approximated  in  past  studies. 

The  jet  intake  can  generally  be  modeled  as  an  open  ended  wave¬ 
guide  with  an  obstacle  (the  engine)  some  distance  down  the  waveguide. 
Figure  1-1  illustrates  the  geometry.  The  problem  can  then,  in 
principle,  be  solved  by  the  generalized  scattering  matrix  technique. 
The  scattering  matrices  of  the  significant  scatterers,  if  known, 
can  be  self-consistently  manipulated  to  produce  the  backscattered 
field.  Let 


u^  represent  the  incident  field, 
u^^  =  the  backscattered  field, 

[SiiJ  =  a  matrix,  representing  (in  some  sense)  the  direct 
^  ^  backscatter  of  the  open  waveguide 
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Si2  =  radiation  characteristics  of  waveguide  modes 

$2^  =  coupling  matrix  of  incident  field  to  waveguide  modes 

S22J  =  reflection  of  waveguide  modes  from  the  open  end 

Sfa  ]  =  reflection  of  waveguide  modes  from  the  obstacle 
«  • 

T2|j^  =  transmission  down  waveguide 

Tb2  =  transmission  back  from  waveguide  to  aperture. 


Note  that  off-diagonal  terms  in  each  transmission  and  reflection 
matrix  will  represent  mode  conversion  from  one  mode  to  another. 


From  Figure  1-1  we  can  easily  show  that: 

bs  r  i  .  <• 
u  —  S^i^u  ^*12^2 

ri  "  ^21^^^  * 

*"2  ^  ^b2  ^b  ^2b  '"1 


•"2  ^b2  ^b  ^2b  ^^21^^  S22r2) 

[^■''’b2^b''’2b^22j''2  *  ''’b2V2b^2l“^ 
[^^b2V2b^  ^  ■  ^22]*'2  ^  ^21*^^ 

*"2  '  [^^b2^b^2b^  ^  ■  ^22]  ^  ^21*^ 
"{^11  *  ^12[^^b2V2b^  ^  ■  ^22] 


(1-1) 

(1-2) 

(1-3) 

(1-4) 

(1-5) 

(1-6) 

(1-7) 

(1-8) 


The  matrices  S-.^  S,2.  ^2^  and  $22  have  been  solved  for  by  the  Wiener- 
Hopf  technique.  These  results  are'^quite  complicated  algebraically. 

The  heart  of  the  problem,  however,  is  calculation  of  the  Wiener-Hopf 
factorization  functions.  Considerable  effort  has  been  given  to  cal¬ 
culating  these  functions,  which  is  presented  in  Appendix  A.  This 
study  has  emphasized  circular  waveguides,  because  much  prior  work 
has  been  done  in  the  area,  and  the  synmetry  of  the  geometry  simplifies 
the  problem  somewhat.  Also,  since  engine  geometries  are  circular, 
non-circular  inlets  require  an  additional  model  of  the  mode  conver¬ 
sion  as  the  energy  travels  down  a  waveguide  of  varying  cross-section. 
Thus  our  matrices  T.^  and  T21.  are  diagonal,  having  only  the  relevant 
phase  delay  for  eacR'^mode.  ^This  study  does  not  address  further 
development  of  this  effect  for  non-uniform  waveguides. 
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Summary  of  Circular  Waveguide  Modes 


;  i 
.1 


Propagation  in  a  circular  waveguide  is  limited  to  a  discrete 
set  of  modes,  which  can  propagate  only  for  waveguide  diameters  larger 
than  a  certain  minimum  (cutoff)  diameter.  The  notation  used  in  this 
study  is  consistent  with  Harrington:  TE]]  refers  to  a  mode  with 
electric  field  transverse  to  the  axis  of  propagation.  In  the  Russian 
literature  (Weinstein)  this  is  referred  to  as  a  magnetic  mode,  since 
it  has  a  z-directed  magnetic  field,  and  all  other  field  components 
can  be  simply  derived  from  it. 

The  number  of  modes  that  can  propagate  in  a  circular  waveguide 
goes  roughly  as  the  square  of  the  diameter.  Table  1-1  lists  the  first 

TABLE  1-1 


SUMMARY  OF  LOW-ORDER  CIRCULAR  WAVEGUIDE 
MODES  AND  CUTOFF  FREQUENCIES 


Mode 

cutoff  ^ 

cutoff  ka 

TEll 

.5861 

1.8412 

TMOl 

.7655 

2.4048 

TE21 

.9722 

3.0542 

TM11 

1.2197 

3.8317 

TEOl 

1.2197 

3.8317 

1.3373 

4.2012 

TM21 

1.6347 

5.1356 

TE41 

1.6926 

5.3176 

TE12 

1.6970 

5.3314 

TM02 

1.7571 

5.5201 

TM31 

2.0309 

6.3802 

TE51 

2.0421 

6.4156 

TE22 

2.1346 

6.7061 

TE02 

2.2331 

7.0156 

TM12 

2.2331 

7.0156 

TE61 

2.3877 

7.5013 

TM41 

2.4154 

7.5883 

TE32 

2.5513 

8.0152 

TM22 

2.6793 

8.4172 

TE13 

2.7172 

8.5363 

TE71 

2.7304 

8.5778 

TM03 

2.7546 

8.6537 

TM51 

2.7921 

8.7715 

TE42 

2.9547 

9.2824 

TE81 

3.0709 

9.6474 

TM32 

3.1070 

9.7610 

TM61 

3.1628 

9.9361 

TE23 

3.1734 

9.9695 
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28  modes  in  order  of  increasing  cutoff  diameter.  The  cutoff  is  cal¬ 
culated  by  noting  the  or  J'(j'  )=0  and  then  ka=j'  „  or 

■'  ^  n  nm  n  nm'  ''nm 

'^nm*  ^  f  ~  ^  X  cutoff  in  terms  of  ka  can  be  simply 

derived.  It  seems  more  intuitive  to  discuss  D/X.  Figure  1-3  presents 
electric  field  line  pictures  for  the  first  few  modes. 

This  study  employs  the  5“^“^  time  convention;  the  choice  is 
necessary  for  consistency  with  the  vast  majority  of  Weiner-Hopf 
literature  (Einarsson  et  al).  The  coordinate  system  chosen  will 
be  standard  spherical  and/or  cylindrical  coordinates,  with  the 
origin  at  the  center  of  the  waveguide  mouth.  See  Figure  1-2.  The 
polar  angle  is  6,  the  azimuthal  angle  <|>.  The  radius  will  be  desig¬ 
nated  by  r,  in  spherical  coordinates,  and  p  in  cylindrical  coordin¬ 
ates.  The  waveguide  radius  is  a,  and  the  waveguide  extends  from 
z=0  to  z=-'“  at  p=a. 

Since  a  large  part  of  the  literature  is  concerned  with  the 
Weiner-Hopf  solution  to  scattering  by  a  semi-infinite  circular 
cylinder,  an  heuristic  description  will  be  presented  here.  Tutorial 
discussions  of  the  Weiner-Hopf  technique  can  be  found  in  [Noble  (1958)] 
and  [Morse  and  Feshback  (1953)]. 

The  Weiner-Hopf  technique  is  based  on  the  fact  that  the  fourier 
transform  of  a  causal  function  is  entire  in  a  half  plane.  For  example, 
the  function 

x(t)  t  >  0 

|o  t  <  0  (1-9) 

has  the  Fourier  transform 


X(a))  =  7  e^“^  e"“^  u(t)  dt  =  /  dt 

-00  0 

[e'““’-e°]  =-^  (1-10) 


[Recall  that  we  are  using  time  dependence].  Hence  X(a))  has 

a  single  pole  at  u  =  -  ia  and  is  entire  (has  no  poles)  for  lm(a))>0>-ia. 
In  the  Weiner-Hopf  technique,  the  fourier  transform  of  the  field 
component  is  taken  with  respect  to  the  coordinate  along  one  axis 
of  the  problem  (the  axis  parallel  to  the  semi-infinite  object). 

The  incident  and  scattered  field  are  then  related  by  means  of  the 
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Figure  1-2.  Coordinate  system  for  modeling  engine  inlet. 
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TE  2,1  TE  2,2  TE  2,3 

Figure  1-3.  Electric  field  lines  for  TE  modes 
in  circular  waveguide. 
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Figure  1-3  (continued).  Electric  field  lines  for 
TM  modes  in  circular  waveguide. 


physics  of  the  problem.  Noble  cites  three  different  techniques 
for  doing  this:  field  matching.  Green's  function  integral  equation, 
and  dual  integral  equations.  The  resultant  equation  is  manipulated 
so  that  the  left-hand  side  must  be  entire  in  one  half-plane  and  the 
right-hand  side  in  the  other  half-plane,  with  a  region  of  overlap. 
These  are  then  equated  to  a  third  function,  which  must  be  entire 
over  the  whole  a-plane.  Because  of  Liouville's  Theorem,  this  new 
function  ends  up  being  a  constant,  equal  to  zero.  This  separately 
sets  each  side  of  the  equation  to  zero,  which  leads  to  the  solution 
(by  taking  the  inverse  fourier  transform). 

For  the  cylinder  scattering  problem,  the  incident  field  is 
decomposed  into  cylinder  waves  by  the  Bessel  function  addition 
theorem.  The  transform  is  taken  with  respect  to  the  z-axis: 


e(a,r,(t))  =  f  E(z,r,<p)e'''^^  dz  .  (1-11) 

•oo 


For  each  cylinder  wave,  the  incident  and  scattered  fields  are  related 
by  imposing  the  boundary  conditions  that  the  tangential  electric 
field  must  vanish  at  the  cylinder  walls,  and  the  surface  current 
must  vanish  in  free  space.  These  are  used  to  manipulate  the  problem 
into  a  Weiner-Hopf  form.  This  process  is  described  by  [Weinstein 
(19c9)  and  Einarsson  et  al  (1966)]. 

The  analysis  of  interactions  at  the  mouth  of  an  open  waveguide 
goes  back  to  [Chu  (1940)].  By  using  the  Kirchoff  approximation  the 
radiation  fields  for  the  lowest  order  propagating  modes  for  a  cir¬ 
cular,  and  a  rectangular  waveguide  are  calculated.  In  terms  of  the 
GSMT,  this  would  enable  us  to  compute  elements  of  [Si^]- 

[Levine  and  Schwinger  ( 1948)]  evaluate  the  acoustic  radiation 
and  reflection  characteristics  of  a  hollow  circular  pipe.  An 
integral  equation  for  the  velocity  potential  is  solved  via  the  Weiner- 
Hopf  technique.  Their  study  is  confined  to  symmetrical  modes  incident 
on  the  open  end  of  the  pipe.  In  a  circular  acoustic  waveguide, 
the  symmetrical  modes  are  given  by 


_ are  the  zeroes  of  J  (x) 

om  0 

JqJj  are  the  zeroes  of  Jq(x) 
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where 


[Pearson  (1953)]  was  the  first  to  apply  the  Weiner-Hopf  technique 
to  the  electromagnetic  problem.  He  considers  a  transverse  magnetic 
plane  wave  incident  on  the  open  end,  and  obtains  equations  for  the 
Laplace  transforms  of  the  axial  and  azimuthal  currents  on  the  waveguide 
walls.  These  equations  are  then  solved  by  the  Weiner-Hopf  techniques, 
and  the  currents  can  be  obtained  by  inverse  Laplace  transform.  The 
fields  far  from  the  mouth  down  the  pipe  are  then  evaluated  asymp¬ 
totically  and  found  to  result  from  propagating  modes,  thus  giving 
the  coupling  coefficients.  The  backscattered  fields  are  not  evalu¬ 
ated,  nor  is  the  behavior  near  the  mouth  studied, 

[Jones  (1955)]  analyzes  the  scattering  of  sound  waves  by  a 
solid  semi-infinite  cylinder.  He  considers  both  hard  and  soft  boundary 
conditions  although  numerical  results  are  evaluated  for  only  the 
hard  (du/dn=0)  case.  An  approximation  for  high  frequency  (large 
diameter)  is  developed,  and  a  low  frequency  expression  is  presented, 

A  variational  expression  is  developed  to  establish  limits  of  error 
on  the  approximations  used  in  evaluating  the  exact  expressions. 

Both  the  pressure  field  on  the  cylinder,  and  the  scattered  far  field 
are  evaluated.  The  end  cap  pressure  time  response  due  to  an  incident 
unit  step  is  evaluated  by  taking  the  inverse  Laplace  transform  of 
the  frequency  domain  response. 

[Noble  (1958)]  treats  the  scalar  problem  for  both  radiation  of 
the  lowest  order  mode,  and  coupling  of  an  incident  plane  wave  to 
the  lowest  order  mode  in  a  circular  waveguide.  His  discussion  is 
tutorial  in  nature,  drawing  somewhat  from  Jones'  work.  There  is 
a  more  complete  treatment  of  the  general  technique  used  to  solve 
Weiner-Hopf  problems  than  is  found  in  most  other  references. 

[Jones  (1964)]  applies  the  Weiner-Hopf  technique  to  radiation 
from  a  semi-infinite  hollow  pipe.  This  study  is  somewhat  tutorial 
in  nature,  being  part  of  a  textbook,  and  considers  only  a  few  of 
the  lowest  order  modes.  Jones  notes  that  TE  modes  radiate  more 
efficiently  than  TM  modes.  TM  modes  have  higher  reflection  coef¬ 
ficients  at  the  open  end  of  the  pipe. 

[Einarsson  et  al  (1966)]  give  an  exhaustive  study  of  diffraction 
by  both  the  infinite  and  semi-infinite  circular  cylinder.  The  back- 
scatter  for  a  plane  electromagnetic  wave  incident  on  a  solid,  semi¬ 
infinite,  perfectly  conducting  rod  is  given.  The  scalar  (sound) 
scattering  from  a  semi-infinite  (both  solid  and  thin  walled)  tube 
is  evaluated.  The  radiation  and  reflection  for  a  propagating  scalar 
wave  incident  on  an  open  end  are  evaluated  (reproducing  the  results 
of  Levine  and  Schuinger  and  also  Weinstein).  There  is  a  brief  dis¬ 
cussion  of  finite  cylindrical  resonators  with  one  end  open,  one  end 
closed  (rigid,  Dirchlet  boundary  condition). 
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The  general  solution  for  scattering  of  a  plane  electromagnetic 
wave  from  a  semi-infinite  thin  walled,  perfectly  conducting  tube 
occupies  about  half  the  report.  This  includes  both  the  backscatter 
far-field  [S]]]  and  coupling  coefficients,  for  the  general  plane 
wave  (neither  TE  nor  TM).  The  special  case  of  axial  incidence  is 
considered.  Radiation  from  a  source  inside  the  tube  [S21]  and  re¬ 
flection  from  an  open  end  sre  evaluated,  largely  copied  from 

Weinstein.  Asymptotic  radiation  of  the  far-fields  is  compared  to 
the  Kirchoff  approximation. 

A  substantial  part  of  this  study  concerns  evaulation  of  the 
Weiner-Hopf  factorization  functions.  A  number  of  forms  of  integral 
expressions  which  exactly  define  them  are  developed.  Power  series 
approximations  to  the  low  frequency  are  developed.  Unfortunately 
these  expressions  are  still  quite  complex. 

In  addition,  only  the  magnitude  of  the  functions  is  approxi¬ 
mated;  the  phase  is  not.  A  large  effort  to  develop  high  frequency 
approximation  yields  some  useful  simplifications,  but  the  results 
are  expressed  in  terms  of  another  unknown  function,  albeit  much 
simpler.  Numerical  data  is  given  for  varying  values  of  ka  with  the  a 
parameter  fixed.  This, is  inconvenient  since,  generally,  a=kcose, 
and  we  are  primarily  interested  in  fixing  ka  and  varying  a.  Experi¬ 
mental  data  are  presented  for  finite  cylinders.  No  experimental 
data  are  presented  for  a  semi-infinite  cylinder. 

[Witt  and  Price  (1968)]  analyze  the  problem  of  a  finite  tube 
without  recourse  to  Weiner-Hopf  techniques.  Instead,  the  direct 
backscatter  from  the  rim,  and  the  coupling  coefficients  to  waveguide 
modes  are  calculated  by  what  amounts  to  the  Kirchoff  approximation. 
The  incident  field  tangential  to  the  aperture  is  expanded  as  a  sum 
of  waveguide  modes,  and  a  waveguide  admittance  for  each  mode  is  com¬ 
puted.  The  reradiation  is  also  calculated  via  the  Kirchoff  method, 
using  the  Stratton-Chu  integral.  The  termination  is  modeled  as  an 
impedance,  which  results  in  a  simple  (scalar)  reflection  coefficient. 
The  reflection  coefficient  is  then  transformed  to  its  equivalent 
impedance  as  seen  at  the  mouth  of  the  waveguide.  The  total  waveguide 
admittance  for  that  mode  is  then  calculated,  and  the  scattered  field 
is  expressed  as  a  sum  over  the  waveguide  modes,  times  the  incident 
field  projection  on  that  mode,  times  the  generalized  admittance. 

For  bacjjscatter  with  vertical  polarization  (TE)  (which  corres¬ 
ponds  to  our  (j)  polarization)  they  present  the  formula 


eS(x..o,z')  =  -  M  cose  ^-jkr 


(1-13) 


vJ^(v+Vq) 


I  I  E 


n=0  m=l 


nm 


('')En„(v„)(l+r  ) 

'  nm  0  nm '  V  v  / 

\  'nm  / 


where 

v„  =  kasinS  (incident) 

0  0 

V  =  kasinS  (scattered) 

E  _  =  weighting  coefficient  for  nm-th  mode 

^nm  ”  Pt^opagation  coefficient  for  nm-th  mode 

r _  =  reflection  coefficient  seen  at  aperture  for  nm-th  mode. 

nm 

The  work  of  Weinstein*  in  many  cases  predates  the  works  reviewed 
here.  However,  most  of  Weinsteins's  papers  were  published  in  Russian 
journals.  His  book  [Weinstein  (1969)]  contains  virtually  all  of 
his  earlier  work,  and  is  readily  available.  It  should  simply  be 
noted  that  Weinstein's  work  did  frequently  predate  work  in  the  West. 

[Weinstein  (1969)]  has  collected  all  of  his  earlier  work  on  the 
Wiener-Hopf  technique  into  a  single  volume.  The  book  is  tutorial 
in  nature,  beginning  with  the  simplest  problem,  diffraction  and 
radiation  by  a  plane  parallel  plate  waveguide.  Having  developed 
the  basic  Wiener-Hopf  arguments,  he  proceeds  to  analyze  circular 
waveguides,  first  considering  only  symmetrical  modes  (to  eliminate 
azimuthal  dependence  from  the  problem).  Acoustical  problems  are 
then  solved,  followed  by  the  general  problem  for  electromagnetic 
waves  scattered  and  radiated  by  a  semi-infinite  circular  cylinder, 
including  azimuthal  dependence.  Comparisons  with  answers  obtained 
via  the  Kirchoff  method  are  frequent,  showing  those  circumstances 
under  which  the  Kirchoff  method  works  and  those  under  which  it  fails. 
There  is  a  substantial  discussion  of  the  relative  accuracy  of  Huygens 
principle,  compared  with  edge  diffraction;  a  substantial  point  is 
made  that  edge  diffraction  yields  more  reliable  answers. 

[Kao  (1970)]  presents  a  completely  novel  approach  to  scattering 
from  cylinders.  He  determines  the  currents  on  finite  cylinders  by 
using  point  matching  and  then  calculates  radiation  patterns  from 
the  currents.  For  the  semi-infinite  cylinder  [1970b]  he  sets  up 
two  sets  of  points  with  slightly  different  interpoint  spacing.  By 
various  manipulations  of  these,  he  determines  the  magnitude  of  the 
traveling  wave  launched  on  the  cylinder  by  the  incident  plane  wave, 
as  well  as  the  current  in  the  vicinity  of  the  aperture.  However, 
his  analysis  is  confined  to  broadside  incidence  (6=90°). 

[Bowman  (1970)]  develops  ray-optical  diffraction  expressions 
for  the  scattering  from  the  end  (aperture)  of  the  semi-infinite  wave¬ 
guide,  and  compares  these  with  an  asymptotic  approximation  to  the 


*Also  transtated  as  Vajnshteijn,  Wainstein,  Vainshtein. 


exact  Weiner-Hopf  solutions.  For  the  asymptotic  form  of  the  Wiener- 
Hopf  solution,  he  obtains,  for  direct  backscatter  (6=0) 


-  X  Ip  e 


ikr  ,  ^  ^  J  .m  ^-3/2  gi2mka| 

/irk  a  m=  1  J 


whereas  by  ray  optical  methods,  he  obtains 


(1-14) 


- 


^iii/4  “ 

^ —  ! 

Aka  m=l 


^.mgi2mka 

,m-l  „ 

2  m 


(1-15) 


J 

For  bistatic  scattering  with  axial  incidence  he  obtains  from  approxi 
mating  the  Wiener-Hopf  solution  letting  his  6  go  to  it-6  and  his  (]) 
to  2n-<t> 


.  -j  e^’^'”  2a  cos(kasin0-ir/4) (sin*) 

cos(e/2) 


f,  ^  cos^(e/2)  ?  , 

1^  /au 


m  -3/2  ^i2mka 
m  e 


From  the  ray  optical  approximation  he  obtains 


(1-16) 


-ikr  /  o-  n1/2  cos(kasine-it/4) 


1  + 


giir/4 

Aka 


cos  (9/2)  y 
cos9  ^ 


®  .m  „i2mka 


1  e 


m=l  2^*¥ 


(1-17) 


- 1/7 

In  both  cases  it  is  seen  that  the  term  of  0(ka  '  )  is  identical 
up  to  the  first  two  terms  in  the  summation 


vs 


iei2ka  .  i2e4ika  ,.3gi6ka 
1“  ^372  ,3/2 


^gi2ka  i2gi4ka  ^•3gi6ka 


(l-18a) 


(l-18b) 
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He  attributes  this  difference  to  the  ray-optical  approach,  for  which 
he  first  considered  scattering  by  plane  parallel  plates;  multiple 
scattering  was  effected  by  assuming  a  cylinder  wave  was  generated 
from  each  edge  after  scattering.  This  result  was  then  specialized 
to  having  a  single  point  participate  in  the  scattering  on  each  edge 
but  not  modified  to  take  account  of  the  fact  that  a  cylindrical  wave 
is  no  longer  being  emanated  from  each  edge.  Beyond  this  comment. 
Bowman  does  not  analyze  further,  since  the  exact  result  from  the 
Weiner-Hopf  solution  instructs  us  how  to  modify  the  ray-optic  con¬ 
tribution. 

A  comparison  of  Bowman's  results  with  those  of  Witt  and  Price 
is  illuminating.  If  we  restrict  ourselves  to  on-axis  backscatter 
and  disregard  terms  due  to  reflection  from  the  termination,  the 
formula  presented  by  Witt  and  Price  reduces  to 


E 


s 


00 


=  -jka 


a 


(1-19) 


We  observe  that,  besides  the  opposite  time  convention,  there  is  an 
additional  factor  of  jka  present  in  the  form  presented  by  Witt  and 
Price.  They  remark  that  this  is  recognizable  as  the  physical  optics 
approximation  for  scattering  from  a  flat  conducting  disk.  Since 
this  is  not  in  agreement  with  the  high  frequency  behavior  of  the 
Weiner-Hopf  solution,  we  conclude  that  the  disk  is  not  a  good  high 
frequency  model  for  this  problem.  In  fact,  it  will  be  seen  later 
that  the  disk  does  adequately  model  the  on-axis  scattering  of  a 
cylinder  terminated  by  a  perfectly  conducting  plate. 


[Moll  and  Seecamp  (1970)]  present  a  more  realistic  model  of 
the  engine  geometry.  The  approach  of  Witt  and  Price  is  used  to  model 
the  scattering,  coupling,  and  radiation  at  the  duct  inlet.  The  study 
is  confined  to  TE  modes.  The  termination,  however  is  modeled  by 
two  sets  of  blades,  each  as  shown  in  Figure  1-4,  to  simulate  the 
first  stage  of  a  compressor.  The  blades  are  modeled  as  being  planar 
(normal  to  the  z  axis).  The  two  sets  of  blades  had  different  numbers 
of  blades  and  blade  widths  and  assumed  varying  relative  orientations 
(stator  to  rotor).  The  scattering  at  the  termination  is  modeled 
by  a  similar  procedure  to  that  at  the  inlet.  The  backscattered  field 
is  expressed  as  a  sum  of  modes  traveling  toward  the  mouth.  The  total 
tangential  electric  fields  must  match  at  points  where  there  are 
blades.  For  each  incident  mode,  integrals  were  taken  over  the  area 
covered  by  either  set  of  blades,  forcing  the  fields  to  vanish. 

The  equations  thus  obtained  are  used  to  solve  for  the  scattering 
coefficients  for  modes  generated  at  the  termination.  Then,  the 
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radiation  from  each  mode  is  computed  by  matching  the  radiated  field 
plus  the  incident  field  to  the  internal  field  at  the  duct  aperture, 
where  the  internal  field  is  the  sum  of  the  modes  traveling  down  the 
waveguide  plus  those  scattered  back  by  the  termination.  The  RCS 
is  computed  for  various  orientation  of  the  blade  structure  and  various 
relative  orientations  of  the  two  bladed  structures,  giving  a  range 
of  modulation  of  the  RCS  caused  by  the  rotor  motion.  They  used  31 
and  37  blades  on  the  respective  blade  structures.  There  is  a  brief 
discussion  of  a  non-planar  cap  in  the  duct  termination,  but  the  idea 
is  not  developed. 

[Lee  et  al  (1973)]  are  primarily  concerned  with  the  measurement 
errors  made  when  field  strength  is  measured  with  a  sensor  boom. 

The  effect  of  the  presence  of  the  boom  is  analyzed  via  the  Wiener- 
Hopf  technique,  and  the  relative  distortion  thus  introduced  is  cal¬ 
culated.  The  relevant  part  of  this  paper  deals  with  the  calculation 
of  the  Wiener-Hopf  factorization  functions.  Although  the  general 
factorization  expression  is  developed  (as  an  infinite  product  of 
factors),  the  general  expression  is  clearly  too  complex  to  be  useful. 

A  low  frequency  form  is  developed  for  the  n=l  case  for  both  L+  and 
M+  functions.  The  n=l  case  is  relevant  to  low  frequencies  because 
the  TEii  mode  has  the  lowest  cutoff  frequency,  hence  is  also  the 
slowest  decaying  evanescent  mode  when  all  modes  are  cut  off.  Un¬ 
fortunately  the  low  frequency  expression  presented  do  not  appear 
to  fit  very  well  with  data  computed  in  this  dissertation  by  numerical 
integration  of  the  exact  integral  defining  the  functions. 

(Lee  et  al)  indicate  a  constant  phase  for  the  low  frequency,  but 
it  appears  that  the  phase  varies  rather  rapidly  as  the  frequency 
increases  from  O.C,  The  formulas  do  work  out  to  the  correct  "DC" 
form. 


[Mittra  et  al  (1974)]  present  a  detailed  report  which  includes 
reconciliation  of  the  numerical  solution  with  experimental  data. 

The  complete  Wiener-Hopf  solution  for  the  semi-infinite  cylinder 
is  presented.  There  appear  to  be  some  errors  in  the  results  presented, 
since  the  scattering  coefficients  for  direct  backscatter  (Equation 
(3-86))  could  not  be  reconciled  with  those  presented  by  [Bowman  (1970)]; 
also  the  two  direct  polarization  solutions  for  Sqq  and  Saa  do  not 
reduce  to  the  same  answer  on  axis,  nor  do  they  reduce  to  the  answer 
given  by  [Einarsson  et  el  (1966)]  for  on  axis  backscatter  (Equation 
(5-63)).  The  solution  for  large  pipes  is  simplified  by  an  asymp¬ 
totic  approximation  for  the  L+  and  M+  functions;  these  are  expressed 
as  functions  of  the  series 

I  m-3/2  ei'"(2ka-n7t+it/2)  (1.2o) 

m=l 


and  the  same  summation,  including  either  even  terms  or  odd  terms. 


The  Generalized  Scattering  Matrix  Technique  (GSMT)  is  used 
to  formulate  the  problem.  Considerable  effort  is  expended  expressing 
the  individual  elements  of  the  matrix  in  terms  of  the  Wiener-Hoof 
solution. 

Solutions  for  other  geometries  are  developed  by  the  ray-optical 
method.  The  geomet'^ies  considered  are  an  ellipsoid,  elliptical  plate, 
and  semi-infinite  elliptical  cylinder.  The  semi-infinite  elliptical 
cylinder  is  analyzed  in  terms  of  diffraction  by  the  edge  of  the 
cyl inder. 

These  solutions  are  then  combined  to  estimate  the  total  RCS 
for  an  aircraft.  Numerical  results  were  computed  for  selected  fixed 
frequencies  for  two  aircraft.  For  calculations  involving  a  termination 
in  the  jet  intake,  this  was  modeled  earlier  as  a  perfect  conductor, 
or  as  a  dielectric  plug  of  infinite  depth,  or  as  a  dielectric  plug 
of  finite  depth. 

[Chuang  et  al  (1975)]  present  an  extension  of  the  report  by 
[Mittra  et  al  (1973)].  Starting  with  the  exact  solution  for  bistatic 
scattering  from  a  semi-infinite  pipe,  they  develop  approximate  high 
frequency  expressions  for  the  factorization  function,  based  on  an 
asymptotic  evaluation  of  the  integral  defining  these  functions. 

Both  factorization  functions  can  be  rather  simply  expressed  in  terms 
of  a  third  function,  called  a  modified  Lerch  function  of  order  3/2, 

L(x.v)  .  I  x.-'  (,.21) 

m=l 


They  then  derive  via  the  Mellin  transformation  a  twelve-term  series 
representation  with  complex  coefficients,  which  enables  the  compu¬ 
tation  of  this  function  for  v  =  3/2  to  be  mechanized  trivially. 

The  results  thus  obtained  show  excellent  agreement  with  results  ob-‘ 
tained  by  direct  numerical  integration  of  the  defining  integral, 
down  to  the  cutoff  frequency  of  the  lowest  mode  of  that  order.  Finally 
there  is  a  simplification  of  the  infinite  sum  which  is  present  in 
the  scattering  calculation  mode  by  use  of  the  asymptotic  form  with 
one  of  the  Bessel  function  addition  theorems. 

There  appear  to  be  a  few  errors  in  this  paper.  These  are  summar¬ 
ized  in  Appendix  B. 

[James  and  Greene  (1978)]  indicate  that  both  theoretical  and 
experimental  results  show  substantial  sensitivity  to  wall  thickness. 
They  indicate  that  the  exact  Wiener-Hopf  solution  based  on  infinitely 
thin  walls  breaks  down  when  the  wall  thickness  is  larger  than  .OU. 

The  results  of  their  summary  seems  to  be  that  the  radiation  patterns 
of  thick  walled  pipes  are  narrower  than  given  by  the  Wiener-Hopf 
solution.  Beyond  noting  this  effect,  we  will  not  further  discuss 
this  problem,  since  it  would  require  an  entirely  separate  study. 
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CHAPTER  II 

DISCUSSION  OF  WIENER-HOPF  SOLUTION 


The  exact  Wiener-Hopf  solution  to  the  diffraction  by  a  semi¬ 
infinite  circular  waveguide  can  be  found  in  [Einarsson  et  al  (1966)]. 
The  significant  formulas  from  this  report  are  summarized  in  Appendix 
B.  The  major  numerical  problem  —  that  of  calculating  the  factori¬ 
zation  functions  --  is  discussed  at  length  in  Appendix  A.  This  chapter 
discusses  the  physical  significance  of  some  of  these  formulas.  We 
reiterate  that  the  time  dependence  is  assumed  and  suppressed. 

The  coordinate  system  is  a  substantial  stumbling  block  since 
the  incident  and  scattered  field  are  defined  with  respect  to  §  and  $ 
unit  vectors  which  are  themselves  functions  of  angle.  Further  com¬ 
plicating  matters,  [Einarsson  et  al  (1966)]  and  [Mittra  et  al  (1974)] 
use  different  coordinate  systems.  This  study  uses  the  coordinate 
system  illustrated  in  Figure  2-1,  which  is  the  same  as  [Mittra  et 
al  (1974)],  since  it  forms  a  self-consistant  reference  for  scattering 
calculations. 

The  incident  field  is  assumed  to  come  from  the  angle 
Looking  in  along  the  y-axis  toward  the  origin,  we  see  Figure’2-2. 
Looking  in  along  the  x-axis  toward  the  origin  we  see  Figure  2-3. 

Hence  we  obtain  the  following  unit  vectors. 


§  =  -  X  cos0^.  -  z  sin0^. 


(2-la) 


$  =  -  y  .  (2-lb) 

For  th?  scattered  field,  we  are  not  necessarily  constrained  to  <>5='rr 
(for  bistatic  scattering).  Hence,  we  obtain  general  formulas  for 
the  unit  vectors. 

0  =  X  cosOg  coS(|)g  +  y  cos0^  ^  ^inO^  (2-2a) 


$  =  -  X  sin(j)g  +  y  cos^)^ 


(2-2b) 


In  our  coordinate  system,  the  incident  and  scattered  fields  are  in 
the  same  coordinate  system.  In  [Einarsson  et  al  (1966)],  the  direction 
of  the  z-axis  is  reversed,  but  0^  and  the  unit  vectors  for  the  inci- 
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dent  field  are  defined  in  such  a  way  that  they  physically  point  in 
the  same  direction  as  in  our  system.  Thus  cos0.,  sine,.,  eL  e]  are 
all  the  same  in  Einarsson's  and  our  coordinate  systems.  (See  Figure 


2-2  in  Einarsson;  E^  occurs  for  g=0,  e]  for  g  =  ^  with  a  being  re- 
0  ^  ^ 

placed  by  0.).  In  the  scattered  coordinate  systems,  §  and  $  are 
reversed  in 'sign  (hence  also  Eg  and  E*),  and  cose  has  the  opposite 
sign,  but  sine  is  unchanged,  '’in  addition  (j)  is  replaced  by  -(|) 
since  the  x-axTs  is  the  same  in  the  two  systems. 


A.  On-axis  Results 


The  most  elementary  case  of  the  Wiener-Hopf  solution  can  be 
found  by  taking  the  limit  as  0^0.  In  taking  this  limit  we  will 
consider  both  the  backscatter  and  coupling  of  energy  into  the  wave¬ 
guide. 


Taking  the  limit  of  Equation  (B-3)  (which  gives  the  general 
bi static  scattering  from  the  rim  for  6  incident  and  scattered 
fields)  with  4>5=Tr,  0.=0  =6  (in  the  limit  6->0)  we  obtain 


S 


00 


where 


lim  ^ 


=  Neumann  epsilon  =  l;n=0 
=  2;n?«0 

JnO  ''s  cylindrical  bessel  function 
n  P 

k  is  the  wavenumber  = 

a  is  the  pipe  radius. 


of  the  first  kind 


L+,  M+  are  the  factorization  functions;  these  are  functions 
of  (n,a,k);  k  is  always  understood  and  n  is  usually 
obvious  from  the  context.  Hence  the  only  argument  that 
is  explicitly  given  is  a.  L+(k)  means  L+(n,a=k,k) 


f^  =  nL+(k)/2kaM+(k) 


1  im 
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2(-1)(ka)^ 
4[L+(k)]  ^ 


(2-4a) 
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(2-4b) 


ia^k 

-5-5 - 5 - 5-  (n=l  understood) 

4k^aV(k)^-L+(k)‘^ 


(2-4c) 


since  is  identically  0. 

When  combined  with  Equation  (B-2),  we  obtain 


j'kr 


-ia^k 


4k^aV{k)^-L+(k) 


2  (n=l  understood)  (2-5) 


This  can  be  seen  to  be  the  same  as  Equation  (5-163)  of  [Einarsson 
(1966)]  with  r-^z;  the  same  procedure  applied  to  Sm  (which  gives 
the  general  bistatic  scattering  from  the  rim  for  |  incident  and  scat¬ 
tered  fields)  leads  to  exactly  the  same  result. 
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1-1 


(2-6b) 


4kM+(k)^(l-f^) 


(2-6c) 


-ia^k 


4k^a^M+(k)^-L+(k)^ 


(n=l  understood) 


(2-6d) 


We  see  immediately  from  Equations  (B-5)  and  (B-6)  (which  give 
the  cross-polarized  backscatter)  that  sin  mT=0  so  no  cross-polarized 
field  is  generated. 


The  radar  cross-section  is  defined  by 

,  2 
O 

o  =  lim  4iir 
r+oo 


(2-7) 


Applying  Equation  (2-5),  we  obtain 
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4k^a^M+(k)^-L+(k)^ 


(4k^a^M+(k)^-L+(k)^ 
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4k^aV(k)2-L+{k)^ 


(2-8a) 


(2-8b) 

(2-9) 


The  results  of  Equation  (2-9)  are  evaluated  and  plotted  in  Figure 
2-4.  It  is  of  interest  to  note  that  the  normalized  cross-section 
peaks  at  D/X=.55,  with  o/iTa'^=10.4  dB,  and  then  decreases  very  rapidly. 
The  lowest  order  propagating  mode  (TEll)  is  enabled  at  D/X=.586; 
hence  this  peak  occurs  just  below  cutoff  for  this  mode,  and  as  the 
mode  is  able  to  transport  energy,  the  cross  section  rapidly  decreases. 
The  first  three  peaks  and  corresponding  modes  are  summarized  in  Table 


2-1. 


TABLE  2-1 

SUMMARY  OF  ON-AXIS  RCS  PEAKS  AND  MODE  ACTIVITY 


Location  of 
peak  (D/x) 

Height  of 
peak  (dB) 

Cutoff  for 
mode  (D/x) 

Mode 

.55 

10.4 

.586 

TEll 

1.68 

3.5 

1.697 

TE12 

2.69 

■  2.5 

2.717 

TE13 

The  rapid  variation  of  cross  section  in  these  regions  is  related 
to  the  strong  coupling  of  the  normally  incident  plane  wave  to  the 
TE-j  modes,  as  discussed  below. 

[Bowman  (1970)]  and  [Chuang  et  al  (1975)]  use  their  asymptotic 
forms  for  L+  and  M+  to  obtain  a  simple  form  for  the  on-axis  cross 
section. 
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This  leads  to 
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(2-11) 
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Equation  (2-11)  is  plotted  in  Figure  2-5  and  yields  good  results 
for  D/X>1. 

Taking  the  inverse  Fourier  transform  of  Equation  (2-5)  with 
the  spatial  dependence  suppressed,  we  obtain  the  time  domain  response, 
shown  in  Figure  2-6.  This  represents  the  backscattered  field  as  a 
time  function  resulting  from  a  normally  incident  plane  wave,  impulsive 
in  time.  This  result  was  generated  via  the  discrete  Fourier  trans¬ 
form  of  the  backscattered  formula  for  0<D/X<3.2.  Naturally  the  time 
domain  response  is  dominated  by  the  peak  at  D/X=.55,  which  results 
in  damped  oscillations.  The  time  domain  response  can  be  related 
(perhaps  dubiously  at  low  frequencies)  to  a  ray  optic  model  of  scat¬ 
tering  at  the  mouth.  The  optic  model  was  developed  by  [Bowman  (1970)]. 
For  example,  if  we  take  the  simplest  asymptotic  approximation  for 
the  L+  and  M+  functions,  given  by  Equations  (A-27),  (A-29),  (A-31), 
(A-58)  and  (A-59),  and  substitute  them  into  Equation  (2-5),  suppres¬ 
sing  the  spatial  decay  and  propagation, 
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E^  ^  1  -4^ia^k _ 

E^  ^  4k^a^M+(k}^-L+(k)^ 


(2-12) 


-1/2 

retaining  only  terms  of  0(1)  and  0(ka  ),  with  a  unit  incident 
field 
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(2-13) 
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giTT/2(m+l/2)g^  c 


(2-14) 


Taking  the  first  few  terms  of  the  summation  explicitly 
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Figure  2-6.  Inverse  Fourier  trans 
backscatter  with  weight' 
Gibbs- type  ri 
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Note  that  e  ^  results  in  a  time  delay  of  ■—  .  We  employ  the  Fourier 
transform  pair 

J_  ^  Mill 
✓w  /2Trt 

and  its  Hilbert  transform 

1.  < 

/-2irt 

Taking  the  transform  of  the  frequency  domain  expression 
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leads  to 
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This  is  plotted  in  Figure  2-7. 

We  note  in  Figure  2-6  that  the  sharp  peaks  at  t=m-|^  become 

almost  negligible  after  m=3,  and  the  response  shown  in  the  discrete 
Fourier  transform  results  seems  to  be  dominated  by  the  ringing  as¬ 
sociated  with  the  lowest  frequency  resonance.  Since  this  resonance 
is  not  well  predicted  by  the  asymptotic  forms  (compare  Figures  2-4 
and  2-5)  it  is  not  surprising  that  the  long-time  or  steady-state 
response  is  not  well  predicted.  This  resonance  appears  to  be  related 
to  the  currents  excited  at  the  mouth,  and  to  some  extent,  associated 
with  an  exterior  natural  resonance  of  the  cylinder,  since  all  the 
waveguide  modes  are  evanescent  in  this  region. 
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-igurp  2-7.  Inverse  Fourier  transform  of  asymptotic  on-axis  backscatter. 


For  purposes  of  comparing  the  relative  importance  of  these 
coupling  coefficients,  the  power  flow  associated  with  each  of  the 
modes  can  be  computed.  For  TE  modes,  with  an  incident  plane  wave 
of  unit  amplitude  based  on  [Collin  {I960),  p.  179]  the  power  flow 
is  computed  as 
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We  assumed  a  unit  incident  field,  which  has  an  incident  power 
1  i  ^  1 

density  of  ^  normalized  power  flow  becomes 

(normalized  to  the  average  power  flow  and  the  area  of  the  waveguide 
mouth) 
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Similarly,  for  TM  modes,  we  obtain 
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The  power  flow  coefficients  are  plotted  in  Figure  2-8  for  the 
first  three  and  modes,  for  axially  incident  plane  waves. 

The  reason  for  the  behavior  of  the  on- ax  is  cross  section  becomes 
somewhat  clearer,  since  we  see  that  the  TE-i  modes  all  couple  much 
more  strongly  to  the  axially  incident  plane'"wave  than  the  TM,  modes. 
Hence,  the  region  in  which  they  are  enabled  exhibits  a  much  more 
dramatic  variation  due  to  the  size  of  the  change  of  power  absorbed. 
[Weinstein  (1969),  p.  151]  discusses  this  behavior  and  proves  that 
asymptotically  the  optical  cross  section  equals  the  sum  of  the  ab¬ 
sorption  cross  sections  for  the  TE,  modes.  Figure  2-8  can  be  com¬ 
pared  with  Weinstein's  Figure  53. 
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Figure  2-8.  Coupling  of  on-axis  incident 

modes  -  normalized  p 


B. 


Off-axis  Behavior 


The  radiation  patterns  of  the  various  waveguide  modes  have 
been  produced  in  the  literature;  all  that  is  produced  here  are  a 
few  curves  for  comparison.  Figures  2-9  and  2-10  with  ka=4  can  be 
compared  with  [Weinstein  (1969)]  Figures  46  and  48.  Figure  2-11 
and  2-12  with  ka=12.77  can  be  compared  with  [Narasimhan  (1979)]. 

Of  more  practical  importance  is  the  relationship  of  these 
radiation  patterns  to  radar  cross-section.  As  is  discussed  in  Section 
II-C,  the  radiation  pattern  of  a  single  waveguide  mode  in  a  certain 
direction  can  be  seen,  by  reciprocity,  to  be  directly  proportional 
to  the  coupling  of  an  incident  plane  wave  from  that  direction  to 
the  appropriate  waveguide  mode.  If,  for  the  moment,  we  assume  that 
the  cylinder  is  terminated  with  a  perfectly  conducting  flat  wall, 
each  propagating  mode  will  be  reflected  back  to  the  waveguide  mouth 
unattenuated.  For  a  monostatic  radar  system,  this  will  mean  that 
at  some  angle,  which  happens  to  couple  well  to  a  particular  waveguide 
mode,  that  waveguide  mode  will  radiate  equally  well  in  that  direction. 
Therefore,  it  becomes  of  considerable  importance  to  know  the  directions 
and  relative  strengths  of  the  coupling  coefficients  of  the  various 
waveguide  modes.  Obviously,  for  6=0,  we  can  see  from  Figure  2-8 
that  the  TEi]  mode  dominates  all  others.  From  Table  1-1  we  see  that, 
for  increasing  waveguide  diameters,  the  TMq^.  TE^],  TM]],  TEqi  are 
the  next  modes  to  propagate.  In  a  manner  analogous  to  Figure  2-8, 
the  relative  importance  to  radar  cross-section  for  these  modes  is 
shown  in  Figure  2-13.  However,  Figure  2-13  differs  in  that  the 
direction  is  varied  for  each  mode  so  that  the  incident  field  is 
assumed  to  come  from  the  optimum  direction.  This  direction  is  plotted 
in  Figure  2-14,  For  example,  for  D/x=1.5,  the  TEqi  has  a  maximum 
for  E(j,  at  30°,  the  TMq]  and  TE21  modes  also  have  maxima  at  30°,  but 
for  E0.  However,  the  relative  magnitudes  are  0.7,  -6.8,  and  1.0 
dB,  respectively.  Figure  2-13  might  be  termed  the  cavity  cross- 
section,  since  it  is  the  contribution  to  the  RCS  which  coupling  to  . 
and  radiation  from  a  single  mode  would  produce,  assuming  that  only 
that  mode  contributed  to  backscatter.  Naturally,  this  is  not  true, 
as  other  modes  will  contribute  to  some  extent,  as  well  as  the  direct 
backscatter  from  the  rim.  We  see  that  for  on-axis  scattering,  the 
"cavity  cross-section"  produced  by  the  TE-.-.  mode  is  remarkably  similar 
to  that  of  the  disk. 

For  an  incident  plane  wave,  we  have  from  Equations  (B-9)  and 
(B-10)  (which  give  the  coupling  of  incident  plane  waves  to  the  axial 
waveguide  field) 


'  V^i 

/F  v^i 

1  *  E  Z  1  ) 

... 

(2-24) 


34 


2-9.  Radiation 


1.273) 


Figure  2-11,  Radiation  pattern  for  TEni  mode  with  ka=12.77  (D/X=4.065 


.77  (D/X=4.065) 


Figure  2-14.  Angle  for  beam  maximum  for  six  lowest-order 

waveguide  modes. 
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From  Equation  (B-17),  (which  give  the  radiation  patterns  for  the 
waveguide  modes)  we  obtain  scattered  fields  of  the  form 


CegOs.t^s) 


Sh^®s’*s^ 


(2-25) 


Consider  first  the  topline  of  Equation  (B-9). 


E^ 


Anm(ei.*i)cos  n,^ 

n'  nm 


-  ■> 

nm 


(2-26) 


A  a 

This  gives  us  E°jjj  =  A^(e^,(ti^)  in  Equation  (B-15)  except  for  the 

direction  of  propagation.  If  we  assume  perfect  reflection,  and 
ignore  the  sign  change,  (since  we  are  interested  only  in  magnitude) 
this  leads  to 


ikr 


(2-27) 


The  RCS  becomes 


0  =  1  im  4'trr 
r-«» 


=  4Tr 


(2-28a) 

(2-28b) 


Since  we  chose  4l^•=^r,  we  likewise  set  (|)^=tt  and  vary  e=e^.=ej  to  its 

optimum  value.  Note  that  for  computations  involving  and  8®^^, 

there  is  a  sin  nij)  dependence  in  the  modal  fields.  To  adjust  this 
to  conform  with  Equations  (B-15)  and  (B-16),  it  is  necessary  to  rotate 
the  coordinate  system.  This  results  in  radiation  patterns  evaluated 

at  4)j  s  Hence,  we  obtain,  for  the  remaining  three  cases: 


0  =  4it 


SH<®s’»s)r 

♦5  =  ti/2,  (|>^  =  ir  ,  Gj  =  0^.  =  e 


(2-29) 
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o  =  47r  C^E(es.4>s)  ^ 

2 

$5  =  '»f/2,  4>i  =  IT, 

0g  =  0.  =  0  (2-30) 

2 

♦  s  ■  -  i>.  85  - 

0.  =  0  (2-31) 

The  backscatter  directly  from  the  rim  is  discussed  by  [Chaung 
et  al  (1975)].  They  used  a  pulse  radar  and  extracted  the  RCS  of  the 
rim  based  on  the  first  return  pulse.  The  agreement  between  theory 
and  experiment  is  better  in  some  cases  than  others,  but  fails  to  match 
the  detailed  pattern.  The  normalized  RCS  can  be  derived  from  Equations 
(B-1)  through  (B-4)  as  follows. 

s  i 

^00  '^0  r 

(2-32) 

2  ^0 

0  =  lim  4iTr  —j 

Eg 

(2-33a) 

=  4^|Se0|^ 

(2-33b) 

0  _  4  |<.  |2  _  |2  ^ 

;;?-^Pee|  ^^001 

(2-34)- 

Similarly 

^  "  1> 

(2-35) 

Comparison  of  Figures  2-15  and  2-16  with  Figures  2  through 

5  of  [Chuang  et  al  (1975)]  reveals  that  despite  use  of  the  complex 
scattering  form,  and  more  accurate  computation  of  the  factorization 
functions,  some  substantial  discrepancies  still  exist  between 
theory  and  experiment.  Possibly  the  explanation  is  that 
any  real  cylinder  must,  of  necessity,  have  a  finite  wall  thickness. 
[James  and  Greene  (1978)]  showed  that  this  leads  to  substantially  dif¬ 
ferent  results  than  obtained  with  the  infinitely  thin  "knife-edge". 

Of  course,  the  accuracy  of  the  measured  results  is  a  possible  source 

of  discrepancy. 
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Figure  2-15.  Rim  backscatter  cross-section  for  semi-infinite  cylinder 

with  ka=7.261  {D/X=2.311). 


Figure  2-16.  Rim  backscatter 


iroci ty  Considerations 


It  is  evident  from  physical  considerations  that  there  must 
be  some  relationship  between  the  coupling  of  an  incident  plane  wave 
to  a  waveguide  mode,  and  the  waveguide  radiation  pattern  of  that 
mode.  This  will  be  shown  by  the  reciprocity  theorem. 

From  [Harrington  (1969),  p.  116-117],  if  the  current  produces 
fields  (E®,H^),  and  current  produces  fields  (E^,H^),  then 


-  X  -  E^  X  H^)*ds  =  -  ?.J®)dv  (2-36) 

s  V 

where  the  surface  and  volume  are  of  finite  extent.  Generally  speaking, 
reciprocity  is  applied  to  sources  and  matter  of  finite  extent.  In 
this  case,  however,  the  matter  is  of  semi-infinite  extent.  Usually, 
it  is  shown  that  the  surface  integral  vanishes  far  from  all  matter 
and  sources.  In  this  case,  we  will  assume  that  it  vanishes  external 
to  the  waveguide,  and  thus  need  only  evaluate  the  'power  flow' 

(E®  X  H^)  inside  the  waveguide. 

First,  we  apply  reciprocity  to  determine  what  source  will  produce 
a  plane  wave  incident  on  the  waveguide.  We  assume  that  there  exists 
a  current  d2pole,  impulsive  in  space,  of  either  0  or  $  orientation, 
located  by  Rg  =  (0=00*  '"=Ro),  where  Rq  »  a,  R  »  X,  as  shown 

in  Figure  2-17.  This  dipole  produces  a  far-field  at  the  origin  of 


(2-37a) 


(2-37b) 


(2-37c) 


for  0  or  $  source  unit  vectors,  respectively.  Hence,  to  produce 
a  'plane  wave'  of  unit  amplitude  in  the  vicinity  of  the  origin,  it 
is  necessary  to  set 


incident  "plane-wave". 
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Let  us  first  consider  TE  modes.  These  have  fields  given  bj 
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(2-38b) 

(2-39) 

(2-40) 

(2-41a) 

(2-41b) 

(2-41C) 

(2-41d) 

(2-41e) 
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+ia^ii,z 

where  e  means  the  fields  propagate  toward  the  waveguide  mouth 

-ia'  z 

e  means  the  fields  propagate  toward  - 

Let  the  b  source  be  located  within  the  guide  at  z  =  -L,  as  shown 
in  Figure  2-18. 


/e“  ioi' 
-  2  \— 

V11„  “JI 


-  na 


n''^nm'  \cos  nij) 


6(z+L)  (2-42) 


This  current  source  will  generate  the  appropriate  TE  mode  of  unit 
amplitude  propagating  toward  the  waveguide  mouth.  It  will  generate 
the  same  mode,  of  equal  amplitude,  propagating  toward  the  infinite 
recesses  of  the  waveguide.  Note  that  since  the  leading  sign  on  H- 
and  H(j,  changes  with  direction  of  propagation,  the  discontinuity  in 
tangential  H  exactly  matches  the  surface  current.  However,  E  E. 
and  are  all  properly  continuous.  ^ 

Next,  according  to  Equation  (B-IO),  an  incident  plane  wave 

with  ()  polarization  produces  a  field  inside  the  guide. 
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’  ^0  n  m  cos  n^ 
nm 
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sin  J„(j'  P/a)  -ia'  z 


(2-43) 
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According  to  Equation  (B-17),  the  TE  mode  propagating  toward  the 
waveguide  mouth  produces  a  far-field  pattern 


Plugging  the  currents  and  fields  into  the  reciprocity  equation, 
taking  the  surface  and  volume  indicated  in  Figure  2-18,  we  obtain 
the  following  result.  For  §  incidence,  h5  has  sin  n(i)  dependence, 
so  we  take  the  bottom  line  of  Equations  (^-41)  and  (2-42). 
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(2-45) 


We  invoke  the  orthogonality  of  modes  to  assert  that  the  sum¬ 
mation  reduces  to  a  single  term,  as  this  integral  vanishes  except 
when  the  mode  and  current  have  the  same  radial  wavenumbers  and  azi¬ 
muthal  dependence.  When  this  is  satisfied. 
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(2-46). 


/;/  f'-T^dv  =  p  •  2 

V  ''nm  ’  "0 


15  iS! 

»  ^0  "^nm 


la'  a  -la'  L 
nm  nm 

e  •  IT 


f  /na 

I  -  'ts’ 


n^nm^ 


pdp  (2-47a) 


2iiB°„kaa, 
nm 


rT?T7  2Vir 


'^n^'^nm^  ^nm 


0  0 


2 


nm 


PI'^P 

(2-47b) 


50 


t  a  » 

—  f  pj  (i'  /a)  dp 


{2-47C) 


2TrB„„kaa'  a 
nm  nm 

TT'  ^ ^ 

''nm 


itB®  ka'  a^ 
nm  nm 

''nm 


Since  is  generated  by  a  field  with  a  sin  dependence, 
rotate  the  radiated  field  generated  by  C.m  by  tt/2.  Hence, 
ate  Cg^  with  a  cos  n4)  dependence.  ^ 


dv 


y- 


This  surface  integral  within  the  waveguide  becomes 
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We  observe  that  H®  = 
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Then,  since  each  of  E^.  Ea,  Hp,  and  are  based  on  H  ,  the  fact 

iples  forces^the  integral  to^vanish  identi- 


that  they  are  scalar  Sulttp 
cally.  Thus  we  are  left  with 
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which  can  be  seen  to  be  satisfied  by  some  simple  algebraic  manipu¬ 
lations. 


For  a  $  incident  plane  wave,  we  note  that  h|  has  a  cos  n(j)  depen¬ 
dence,  so  we  take  the  top  lines  of  Equations  (2-41)  and  (2-42)  for 
both  the  current  and  waveguide  fields.  This  results  in  exactly  the 


same  integrals  for  both  E^-J^  and  with  e-v^,  so  that 
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which  is  satisfied  identically. 
For  TM  modes,  we  begin  with 
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6(z+L)  (2-53) 


This  will  generate  the  appropriate  field  traveling  toward  the  wave¬ 
guide  mouth.  However,  to  generate  the  discontinuity  in  tangential 
H,  the  leading  sign  on  all  terms  must  be  reversed  for  the  mode 
traveling  toward  the  infinite  recesses  of  the  tube.  However,  as 
before,  this  is  of  little  importance  since  it  contributes  nothing 
to  the  reciprocity  integral. 

The  incident  plane  wave  generates  fields 
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For  the  0  incident  field  we  obtain 


o  .0  •  ™ 

—  -r  0  2Tr  a  A„_ia__ae 

rrr  t-dib^  r  r  r  nm  HITI 

JSI  E  .J  d=  /  J  J  - .  J.v,--)  ■ 

V  _oo  0  0  ''ntn  n^'^nm' 


-pj;,{jnmP/a)cos  nd) 


+  $  pa - ^LJE -  sin 

•^nm  p 


ika  e 


■ia__L 

nin 


0  ''ntn  n"'nm' 


t  na  Jn'jnmf/")  ... 

(p  I  -  S  in 

'^nm  ^ 


in  nd)  -  p  6(z+L) 


,  S  .6  '“P 

V  ^pff  J 


pdpdtjtdz 


,  ia__L  I 

nm!  e  ika  g‘  ntn*- 

■^nm  '^nm 

n '  nm 


2  .  n^a^ 


(2-56a) 


0  \  Jnm 


(2-56b) 


_  B  2 

^Vr  C-  ^  ■’n'Jn.P/")'  -to 

^  0  Jr.™  Jl(jnm)  0 

nm  n '  nm ' 


(2-56c] 


fl  ka„„a 

2tt/-^  a® 
v  u .  nm  J  2 1 


ii/j  at  z  n"'nm' 


(2-56d] 


(2-56e) 


ij 


-T  — r  ^ikr  /eT  i4TiR„  -ikR- 

///  ?. J^dv  = ;//  (ecg^^c^E)  V  •  Vu^  -nr^  ^  °6(?-iT^)dv 

{2-57a) 
(2-57b) 


/  ”0  r  4ff1 

■•'  ^GE  nr 


Equating  these  leads  to 
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Finally,  for  4)  incidence,  the  field  generated  has  sin  ncj)  depen¬ 
dence,  forcing  us  to  the  bottom  line  of  the  current  and  field  distri¬ 
butions  (Equations  (2-52)  and  (2-53)).  The  integrals  all  end  up  the 
same,  resulting  in 
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Therefore,  we  have  demonstrated  the  relationship  between  coupling 
of  incident  plane  waves  to  waveguide  modes  and  the  radiation  patterns 
of  waveguide  modes. 
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CHAPTER  III 

SCATTERING  BY  A  REPRESENTATIVE  ENGINE-LIKE 
OBSTACLE  IN  A  WAVEGUIDE 


This  chapter  is  concerned  with  scattering  by  a  simplified  model 
representing  an  engine  geometry  in  the  intake  duct.  The  model  consists 
of  an  infinite  circular  waveguide  housing  an  axially  symmetric  cone 
centered  on  a  flat  plate,  both  perfectly  conducting,  as  shown  in 
Figure  3-1.  This  model  is  based  to  some  extent  on  observation  of 
the  Pratt-Whitney  J-57  turbofan  jet  engine,  on  display  at  the  Air 
Force  Museum,  Wright-Patterson  AFB,  Ohio.  In  this  study,  it  is 
assumed  that  the  TE^i  mode  is  incident;  similar  procedures  could 
be  followed  for  any  other  incident  mode.  The  coordinate  system  is 
centered  at  the  base  of  the  cone,  as  shown  in  Figure  3-2. 

The  calculation  of  fields  scattered  by  the  cone  in  situ  would 
be  known  completely  if  one  knew  the  exact  currents  on  the  obstacle. 

The  scattered  fields  can  be  obtained  straightforwardly  from  the  dyadic 
electric  Green's  function  and  the  dyadic  magnetic  Green's  function, 
based  on  the  integrals  given  by  [Tai  (1973),  p.  9]. 


E(R)  =  icopg  ///  ggi(R.R')-'3(l?')  dv'  (3-1) 

Tf(R)  =  ;jj  ■ 


Although  the  source  singular  term  for  the  dyadic  electric  Green's 
function  is  still  being  discussed  in  the  literature  by[Yaghjian  (1980)] 
and  others,  the  Green's  functions  at  points  not  near  the  source  can 
be  obtained  unambiguously  from  residue  series  expansions.  These 
expansions  consist  of  a  double  summation  over  n(azimuthal  index) 
and  m(radial  index)  of  terms  representing  TE  and  TM  fields,  both 
propagating  and  evanescent.  Hence  the  coupling  to  a  given  mode  can 
be  determined  by  a  single  integral,  allowing  one  to  obtain  either 
E^  or  H^,  based  on  known  currents.  Mathematically  this  can  be  ex¬ 
pressed  as  follows:  Omitting  the  source  singular  term,  one  can  write 
(using  Tai's  notation) 
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Figure  3-1.  Basic  geometry  for  engine  scattering  model 


Figure  3-2.  Coordinates  system  and  dimensions  for 
scattering  from  an  axial  conducting  cone. 
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Tr(h)  represents  TM  electric  fields  or  TE  magnetic  fields;  If  contains 
an  axial  component 

W(h)  represents  TE  electric  fields  or  TM  magnetic  fields  so  that 
W  is  purely  transverse. 
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For  example,  for  a  single  TE  mode,  we  can  obtain 
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z,n,m 


le 
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4Tiy^k  I 

p  y 


///Me^y{+k^)-^(R')  dV  (3-5) 


and  similarly  calculate  E  using  5  The  elements  of  G  ■,  and  G  ^ 
are  written  out  in  Appendix  C. 


Because  of  the  difficulty  of  and  restrictions  on  previous 
solutions  for  current,  a  new  technique  of  solving  for  currents  was 
explored.  A  few  cases  have  been  previously  solved  with  great  effort. 

For  example,  [Tesche  (1972)]  solved  the  problem  of  a  skewed  wire  between 
parallel  plates  by  using  images.  [Wang  (1978)]  used  the  dyadic  elec- 
trir.  Green's  function  to  solve  for  the  currents  in  an  arbitrarily 
shaped  dielectric  body  inside  a  rectangular  waveguide.  [Harrington 
(1961),  pp.  402-406]  gives  a  variational  technique  to  find  stationary 
formulas  for  scattering  using  approximate  current  distributions. 

He  applies  this  to  a  post  in  a  parallel-plate  guide.  Unfortunately, 
this  method  requires  computing  the  self-reaction  of  the  assumed 
current  distribution,  which  in  most  cases  is  equally  as  difficult 
as  solving  the  problem  exactly,  since  it  requires  computing  the 
electric  field  generated  by  the  assumed  current  in  the  source  region. 
Hence,  no  computational  advantage  is  obtained  over  solving  directly 
for  the  actual  current  by  the  method  of  moments,  or  some  similar 
technique. 


In  principle,  the  problem  can  be  solved  by  the  method  of  moments 
by  assuming  the  current  distribution  to  be  a  collection  of  a  series- 
of  pulses  with  unknown  weighting  coefficients,  calculating  the  elec¬ 
tric  fields  generated  by  these  pulses,  and  adusting  the  coefficients 
such  that  the  tangential  electric  field  vanishes  as  the  surface  of 
the  obstacle.  The  problem  with  this  approach,  as  pointed  out  by 
is  that  for  coplanar  field  and  source,  the  dyadic 
s  function  converges  extremely  slowly,  if  at  all, 
in  addition  to  the  problems  involved  with  the  source  singular  term. 
Convergence  of  the  residue  series  is  enforced  by  the  axial  propagation 
icL  |z-z' I  -a__|z-z'| 

factor  e  ,  which  becomes  e  ,  for  large  enough  m. 

For  coplanar  source  and  field,  this  factor  becomes  unity.  In  Wang's 
case  the  problem  was  solved  by  summing  the  resultant  series  (without 
the  convergence  factor)  in  closed  form.  In  the  case  of  the  circular 
waveguide,  we  have  a  Fourier-Bessel  series  whose  summation  in  closed 
form  is  not  readily  apparent.  For  example,  the  z  component  of  G  ^ 
is  given  by 


[Wang  (1978)], 
electric  Green' 
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Taking  asymptotic  forms  for  j  and  J^{  ),  assuming  n=l,  for  coplanar 
source  and  field  point,  we  obt^ain 


J]m  ('n+%)TT 

Jl(x)  cos(x-3Tr/4) 
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z'^  £os({)^cos(|)  j  ll!!tkIcos({nH-Js)Trp/a-3Ti/4)  cos(  (m+}s)Tip 7a-3TT/4) . 
*  k  a  m=l  /pp' 


(3-8) 

It  seems  that  even  if  p  and  p'  are  different,  this  series  does  not 
converge.  The  series  does  converge  for  real  objects  after  inte¬ 
grating  over  volume  currents  of  finite  extent,  since  integration 
over  p'  results  in  dividing  the  asymptotic  terms  by  (m+>i)TT,  and 
integration  over  z  multiplies  the  element  by  the  differential  dz, 
forcing  the  coplanar  elements  to  make  an  infinitesimal  contribution. 
The  reminder  of  the  integral  is  well  behaved,  due  to  the  convergence 
factor  for  noncoplanar  points. 

This  extremely  cumbersome  process  seemingly  cannot  be  avoided. 
However,  poor  convergence  in  the  source  region  suggests  a  similar 
problem  in  free-space  scattering,  and  the  superior  convergence  of 
the  magnetic  field  integral  equation  (MFIE)  over  the  electric  field 
integral  equation  (EFIE)  for  most  obstacle  scattering  P'^oblems^^  ^  1 
We  therefore  proceed,  analogous  to  the  MFIE  to  force  Ug  =  n  x 

just  outside  the  surface  of  the  conducting  body  and  n  x 
just  inside  the  surface.  As  with  the  MFIE,  taking  the  mean  at  the 
surface  results  in  dividing  the  current  by  two.  Assume  the  current 
consists  of  a  collection  of  pulses  of  arbitrary  weighting 


^  »kP<''  -  '•k>“k  (3-9) 

whefVe  u.  is  a  unit  vector  tangent  at  the  surface  at  r.  ,  and  p(r-r.  ) 
is  localized  near  r^.  Then  ^ 
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(3-10) 


Enforcing  the  boundary  condition 


=  n  X  =  n  X  H^'  +  n  X  H^ 


n  X  =*5jg  -  n  X 


n  X  H  (F)  =  J  a|^|sLi|^-n  x  // 


(3-11) 


We  can  then  enforce  this  condition  as  many  times  as  necessary  to 
obtain  the  required  number  of  equations  to  solve  for  the  Sj^'s.  The 
advantage  of  this  approach  is  found  only  by  careful  examination  of 
the  convergence  properties  of  G^^  and  Gjjj2*  The  least  convergent 

term  in  Gg-j  goes  asymptotically  as  m  sin  mx,  but  the  least  convergent 
term  in  Gjjj2  goes  as  1  sin  mx.  After  integrating  over  the  surface, 

we  obtain  ^ or  — — ^ ,  which  are  convergent.  Hence,  we  can 

work  with  surface,  rather  than  volume  currents,  and  construct  a 
relatively  simple  code.  Additionally,  there  appears  to  be  no  dispute 
concerning  the  behavior  of  §2  in  the  source  region. 


With  the  geometry  shown  in  Figure  3-2,  we  define 
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Assume  the  current  consists  of  the  sum  of  pulses 
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(1  n 

P(p-Pk)  +  \  $  p(p-p|^) 


(3-13) 


where  p(p)  is  a  triangular  pulse  function.  This  current  distribution 
is  shown  in  Figure  3-3.  Then  using  Equation  (3-11),  we  obtain 


n  X  =  y  a.  t5sp{P-P. )  -  n  x  //  G  .(R,R')-t  p{^?'-I^)  ds’ 


+  I  $*5P(P-P,^)  -  n  X  //  P(R’-'"k) 


(3-14) 

Next,  to  employ  Galerkin's  method,  we  generate  2N  equations  by  multi¬ 
plying  through  by  each  basis  function  and  integrating  over  the  surface 
S,  to  obtain 


u^*n  X  //  H^(f)p(R-r  )  ds  = 
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q  =  1,2, ...N 
i  =  1,2 


^  ^  A 


where  u-j  ®  <j)  U2  =  t. 

Thus  2N  linear  equations  for  the  2N  unknown  coefficients  a^  and  bj^ 
can  be  obtained,  which  can  be  straightforwardly  solved  by  linear  ^ 
algebra. 

In  order  to  further  speed  computation,  elimination  of  unknown 
currents  in  the  backplane  can  be  obtained  by  using  the  method  of 
images.  Instead  of  the  original  problem,  we  introduce  an  image  such 


that  the  image  of  the  cone  has  image  currents  on  it,  and  the  image 
incident  field  is  assumed  to  be  present,  as  shown  in  Figure  3-4. 

This  automatically  forces  nxE^°^^^=0  and  in  the  backplane, 

hence  reducing  the  amount  of  coplanar  integration  necessary.  With 
this  technique,  we  obtain  a  solution  for  the  currents  on  the  cone, 
then  integrate  these  with  5^2  obtain  the  H-field  in  the  backplane, 
hence  the  currents  in  the  backplane.  The  current  distribution  thus 
obtained  is  shown  in  Figure  3-5.  It  would  be  extremely  difficult 
to  measure  this  experimentally,  and  constraints  of  time  and  money 
prohibit  this  verification. 

The  fields  in  the  region  of  the  cone-base  termination  can  be 
approximated  in  the  spirit  of  the  WKB  approximation,  if  the  cone 
is  slender  enough,  and  the  cone  diameter  varies  slowly  enough  with 
axial  position.  At  each  axial  position  on  the  cone,  it  is  assumed 
that  appropriate  TE  fields  exist  for  an  infinite  coaxial  Ij^ne  of 
the  same  inner  diameter.  Since  £2=0,  this  guarantees  n  x  E  =  0  on 
the  cone  and  outer  waveguide  walls.  The  radial  wavenumber  is  then 
computed  by  solving  the  characteristic  equation 


j;(Mb)  Y;(ya) 


Y;(yb)  j;(ya)  =  0 


(3-16) 


for  the  smallest  positive  y,  with  b  the  inner  diameter  and  a  the 
outer  diameter.  The  values  of  y  are  plotted  in  Figure  3-6.  Next, 
the  scalar  wave  equation  for  H,  is  approximated  by  assuming  that 
y(z)  varies  slowly  enough  so  that  we  can  neglect  the  coupling  of 
y(z)  through  the  radial  function.  Suppressing  the  azimuthal  depen¬ 
dence,  let  H  =  R(p,z)  Z(z).  The  scalar  wave  equation 
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is  approximated  by  setting 


=  0 


(3-17) 


R(p,z)  =  0 


(3-18) 


so  that  the  wave  equation  reduces  to 
+  I IR  +  f  2 

io  »  50  J 


Z(2)t(k^-/)R(p)Z(z)+R(p)  -  0. 
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(3-19) 
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Since  it  was  assumed  that  R(p)  was  an  appropriate  combination  of 
Bessel  functions,  this  reduces  to 


4  *  (kV)Z  -  0 


(3-20) 


Since  m(z)  is  a  known  function,  based  on  the  local  radius,  this 
equation  can  be  integrated  numerically  to  obtain  Z(z).  The  results 
of  this  numerical  integration  are  given  in  Figure  3-7.  Finally, 

to  normalize  the  radial  function  properly,  J|R(p,z)|  da  must  be 


constant  locally  at  all  cross-sectional  planes.  Subsequently,  Hp, 
Ha,  and  H  can  be  determined  at  the  inner  conductor,  and  hence  the 
surface  currents  and  J  ,  which  are  shown  in  Figure  3-8. 


There  are  several  1 imitations^to  ^his  approach.  Since  Hp  =  0 
on  the  inner  conductor  but  H  0,  n  •  H  0.  Maxwell's  equations 
are  only  satisfied  approximately,  not  exactly,  since  the  coupling 
of  p(z)  through  the  radial  functions  was  neglected.  The  behavior 
of  the  fields  near  the  tip  of  the  cone  was  approximated  very  loosely 
since  it  was  assumed  that  the  field  consisted  only  of  incident  and 
scattered  TE-j^  modes  immediately  beyond  the  tip. 
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Figure  3-7.  Solution  for  axial  variation  of  fields  for  TE,,  mode 
in  circular  waveguide  with  cone.  L/a=2.0,  b/a=0.5,  ka=i.O. 
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gure  3-8.  Currents  induced  on  cone  in  circular  waveguide  with 
incident  TEu  mode,  L/a=2.0,  b/a=0.5,  ka=2.0,  solution  by 
coaxial  approximation  with  slowly  varying 
center  conductor. 


CHAPTER  IV 

CONCLUSIONS  AND  DISCUSSION 


The  problem  of  predicting  the  radar  cross  section  of  jet  intake 
cavities  in  a  spectral  region  where  the  cavity  aperture  is  of  resonant 
dimension  has  been  approached  via  the  exact  Wiener-Hopf  solution 
for  a  semi-infinite  waveguide.  Exact  (numerical  integration)  RCS 
computations  have  been  presented  for  the  semi-infinite  circular  wave¬ 
guide  for  electrical  diameters  from  zero  to  three  wavelengths.  These 
computations  smoothly  join  high  frequency  asymptotic  results  and 
also  serve  to  define  the  region  of  validity  for  the  simpler  asymp¬ 
totic  models.  Therefore  one  of  the  matrices  in  a  generalized  scat¬ 
tering  matrix  development  of  the  scattering  by  a  finite  loaded  circular 
waveguide  has  been  completed.  In  the  process,  the  nontrivial  relation¬ 
ships  between  two  earlier  studies  of  the  semi-infinite  cylinder  have 
been  developed  and  certain  errors  in  a  publication  from  one  of  these 
studies  have  been  corrected. 

The  coupling  coefficients  relating  an  incident  plane  wave  to 
internal  waveguide  modes  have  been  extracted  from  earlier  studies, 
recast  in  a  more  convenient  form,  and  evaluated  in  the  resonance 
region.  The  predominance  of  the  TE  modes  over  the  TM  modes  in  power 
absorption  for  axial  incidence  has  been  demonstrated.  The  relative 
importance  of  the  five  lowest  order  modes  for  non-axial  incidence 
has  been  demonstrated  graphically.  The  relationship  between  the 
coupling  coefficients  and  the  radiation  patterns  of  the  waveguide 
modes  has  been  established  via  the  Lorentz  reciprocity  theorem. 

The  relationship  thus  demonstrated  is  that  the  radiation  pattern 
and  the  coupling  coefficient  for  a  given  mode  at  a  given  frequency 
and  polarization  are  proportional.  This  is  of  particular  importance 
to  RCS  calculation  since  a  given  waveguide  mode  will  couple  to  and 
radiate  efficiently  in  the  same  direction  and  with  the  same  polari¬ 
zation. 

Modelling  of  the  leading  surface  of  a  jet  engine  as  a  cone 
on  a  flat  plate  inside  a  circular  waveguide,  we  have  computed  the 
currents  on  this  obstacle,  by  using  the  method  of  moments  and  the 
dyadic  magnetic  Green's  function  appropriate  to  the  interior  of  a 
circular  waveguide.  The  reflection  coefficient  for  an  incident  TE-j^ 
mode  has  been  calculated  based  on  these  currents  at  selected  fixed 
frequencies. 

It  therefore  remains  to  apply  these  results  to  Equation 
(1-8)  to  solve  the  radar  scattering  problem  in  a  self-consistent 


fashion.  This  is  done  for  some  illustrative  cases.  First,  fixing 
a  perfectly  conducting  flat  plate  at  ten  radii  down  the  waveguide 
from  the  mouth  as  shown  in  Figure  4-1,  the  RCS  has  been  computed 
while  varying  the  diameter  frcxn  zero  to  1.5  wavelengths.  The  normal¬ 
ized  RCS  is  shown  in  Figure  4-2.  In  Figure  4-3,  the  normalized  RCS 
is  shown  for  the  case  when  the  waveguide  diameter  is  one  wavelength 
and  the  flat  plate  varies  from  two  to  ten  radii  down  the  waveguide. 

Discussion 

The  calculation  of  scattering  from  the  inlet  mouth  and  coupling 
to  internal  modes  has  been  performed  modelling  the  inlet  as  a  circular 
duct  with  infinitely  thin  walls  and  a  knife-edge  rim.  In  practice 
inlet  rims  have  a  finite  wedge  angle,  and  are  very  rarely  circular. 

The  major  inadequacy  encountered  in  calculating  rim  scattering  for 
the  knife  edge  by  ray  optics  proved  to  be  the  underestimate  of  the 
peak  just  above  cutoff  of  the  lowest  order  propagating  mode  (comparing 
Figures  2-4  and  2-5).  Unfortunately,  for  arbitrary  geometries,  there 
seems  to  be  no  simple  way  to  patch  up  ray-optics  to  solve  this  problem 
Since  the  frequency  range  this  occurs  in  is  low,  it  may  be  possible 
to  evaluate  this  region  by  moment  method  techniques.  Although  this 
was  not  done  in  this  study,  we  can  see  that  the  ray  optics  approxi¬ 
mation  does  work  well  provided  we  are  sufficiently  above  cutoff  of 
the  lowest  order  mode.  For  circular  geometries,  this  condition  is 
that  the  diameter  is  greater  than  one  wavelength. 

Because  of  the  reciprocity  relations  demonstrated  in  Section 
II-C,  it  is  clear  that  discussion  of  radiation  patterns  and  coupling 
coefficients  are  redundant.  The  calculation  of  radiation  patterns 
has  been  discussed  by  [Weinstein  (1969),  pp.  139-150].  Weinstein's 
comments  on  the  use  of  Huygen's  princip'e  for  calculation  of  radiation 
patterns  indicates  that  in  the  transition  region  (near  cutoff) 

Huygen's  principle  performs  poorly,  and  worse  for  TM  modes  than  TE 
modes.  As  Figures  2-13  and  2-14  indicate,  the  relative  importance 
of  TM  modes  to  backscattering  from  cavities  is  very  far  below  that 
of  TE  modes,  and  limited  to  angle  far  from  the  forward  direction. 

For  all  practical  purposes,  forward  scattering  in  the  transition 
region  can  be  calculated  using  the  TE]-]  mode,  adding  the  TE2i  and 
TEqi  modes  as  requirements  may  dictate.  The  reason  for  this  can 
be  grasped  physically  by  considering  Figure  1-3.  For  modes  higher 
than  m=l  and  for  all  TM  modes,  the  oscillation  in  sign  across  the 
waveguide  cross-section  forces  the  average  interaction  to  a  very 
small  value,  regardless  of  incidence  angle.  For  m=l  modes,  the  vari¬ 
ation  due  to  cos  n(|)  dependence  can  be  seen  to  lead  to  an  optimum 
angle,  from  which  the  aperture  fields  appear  to  be  nearly  of  the 
same  sign. 

For  non-circular  geometries,  the  primary  conclusion  we  might 
draw  is  that  the  appropriate  low-order  TE  modes  will  dominate  scat¬ 
tering  in  the  forward  direction.  Since  Wiener-Hopf  solutions  are 


Figure  4-3.  Cross  section  of  semi- inf ini te  circular  waveguide 
for  axial  incidence,  D/X=1.0,  with  flat  plate  at  two  to 
ten  radii  down  guide. 


not  available  for  these  arbitrary  geometries,  it  will  probably  be  nec¬ 
essary  to  construct  moment  method  solutions  for  near  cutoff  frequen¬ 
cies.  Higher  frequencies  can  probably  be  adequately  approximated  by 
ray-optic  solutions  as  done  by  [Pathak  and  Huang  (1980)],  for  example. 

The  impulse  response  waveform  given  by  Figure  2-6  indicates 
that  rim  backscatter  is  really  only  significant  for  times  of  t<5a/c. 
Furthermore,  this  short  time  response  can  be  reasonably  approximated 
ray-optical ly,  confirming  our  earlier  statement  that  high-frequency 
scattering  can  be  adequately  modelled  ray-optical ly.  The  low- 
frequency  scattering,  which  is  not  ray-optic,  corresponds  to  the 
long-time  response,  which  is  typically  of  little  interest. 

Having  enumerated  the  difficulties  in  calculating  scattering 
coefficients  for  object  inside  waveguides,  it  becomes  clear  that 
much  work  remains  to  be  done.  The  use  of  the  magnetic  field  dyadic 
Green's  function  and  appropriate  boundary  conditions  provide  some 
improvement  over  use  of  the  dyadic  electric  Green's  function,  but 
it  is  still  very  expensive  computationally.  The  possibility  of 
azimuthal  asymmetry  (for  example,  a  blade  structure)  has  not  yet 
been  addressed.  One  might  model  the  blade  structure  axisymmetrical ly 
by  using  boundary  conditions  in  the  backplane  Ep  =  0  and  J,j)  =  Hp  =  0 
and  at  the  outer  wall  Jp(P=a)  *  -  H,|)(p=a)  =  0.  ^This  allows  radial 
currents  and  azimuthal  E-fields  to  exist,  but  not  radial  E-fields 
or  azimuthal  currents.  It  further  forces  the  axial  current  to  vanish 
at  the  tip  of  the  blades.  For  a  structure  with  28  blades,  for  example, 
it  would  require  an  extremely  large  diameter  waveguide  before  signifi¬ 
cant  azimuthal  currents  could  flow  on  the  blades.  One  problem  in 
implementing  this  study  would  be  that  the  image  procedure  used  in 
Chapter  III  to  eliminate  coplanar  integrations  could  not  be  used, 
and  it  would  be  necessary  to  integrate  coplanar  source  and  field 
points.  Another  approach  would  be  the  use  of  the  free-space  Green's 
function,  to  analyze  the  structure.  In  this  case  it  would  be 
necessary  to  also  set  up  currents  and  enforce  boundary  conditions 
at  the  waveguide  walls,  thus  vastly  increasing  the  number  of  current 
elements  required.  Another  possible  approach  is  to  segment  the 
obstacle  structure  and  use  waveguide  modes  appropriate  to  coaxial 
waveguides  near  the  obstacle.  However,  the  convergence  properties 
of  this  approach  are  unknown. 
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APPENDIX  A 

EVALUATION  OF  WIENER-HOPF  FACTORIZATION  FUNCTIONS 


The  factori2^Jon  functions  L^(a),  and  are  defined  as 

being  functions  which  are  analytic  in  the  upper  half  of  the  complex 
a-plane,  which  also  satisfy  the  equations 

L+(a)L^(-a)  =  iri  ^  (ya)  (A-1) 

M^(a)M^(-a)  =  iri  j;(Ya)Hj,^  ^  •  (ya)  (A-2) 


where  y  -  A2  -  .2  (Ei  narsson  et  al.  Equations  (5-33),  (5-34)). 

The  exact  function  can  be  evaluated  by  numerical  integration  [Mittra 
and  Lee,  (1971),  p.  107,  Eqs.  (9.8)].  The  approximate  evaluation 
of  these  functions  proceeds  in  two  different  regions,  for  which  the 
Bessel  functions  are  approximated  differently.  Section  A-1  deals 
with  exact  numerical  integration.  Section  A-2  deals  with  low-frequency 
approximations.  Section  A-3  deals  with  high-frequency  approximations. 

A-1 .  Exact  Numerical  Integration 


The  exact  expression  for  L_^(a),  is  given  by 

I  1  <*>+ic  jinfiriJ  (ayk^-z^)H^^^(a/k^-z^)]dz 

-  exp|  ^  ^  " 

'  -OOUC 


I 


z-a 


(A-3) 


where  -Im(k)  <  c  <  Im(a)  <  Im(k). 


Before  proceeding  further,  however,  we  note  that  it  would  be  very 
beneficial  to  not  have  to  integrate  out  to  <».  We  note  that,  for 
Re(z)>Re(k),  the  arguments  of  the  Bessel  functions  become  essentially 
imaginary.  From  Abramowitz,  Equations  (9.6.3)  and  (9.6.4) 


■imi 


.  f,  e  ‘  K„(a 


J?: 


(A-4) 
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F/e  21/s 


inn 


(A-5) 


where  K^(  )  and  I^(  )  are  modified  Bessel  Functions. 


Taking  as^ptotic  forms  (Abramowitz  Equations  (9.7.1)  and  (9.7.2)) 


intr  I -  iniT  _ 

Tri  ^  e  '^K^(al4^-k^)  e~^ 


2  e 


(A-6) 


We  are  integrating  In(a^^-k^)  from  k  to  «,  but  would  prefer  to  inte^ 
grate  ln(l)  from  k  to  since  it  is  trivial.  Hence,  we  form  the 
equation 


-ia 


1-+ («)!-+( -a) 


inr3Jp(Y3)Hp^  ^(ya) 


(A-7a) 


✓-ia(k+a)L^(a)  /-ia(k-a)L^(-a)  =  itYaJ^(Ya)Hj|^ ^ya)  (A-7b) 


We  can  evaluate  this  by  the  same  integral,  except  that  the  argument 
of  the  logarithm  becomes  very  close  to  1  for  Re(z)>Re(k).  Taking 
c>0 


/-ia(k+a)L^(a)=expj^  / 


®  An 


Tia\/k^-z^J^  ( a/k^-z^H^  ^  ^  ( a\/k^-z^) 


z-a 


dz 


V  / 

(A-8) 
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Similarly 


/-ia(k+a) ( a) =exp 


1 

2iT 


In 


7ravk^-2^J 


•(a/k^-2^)H, 


2-a 


>dz 


\  j 

(A-9) 


These  integrals  have  the  advantage  that  we  can  essentially  integrate 
on  the  range  where  -(ka+C)  <  z  <  (ka+C)  and  the  total  error  thus 
introduced  can  be  made  negligible  by  proper  choice  of  C. 

The  complex  value  of  k  requires  some  thought.  The  argument 
introducing  this  imaginary  part  is  given  by  [Einarsson  (1966), 
p.  147].  The  correct  answer  is  defined  as  the  one  obtained  in  the 
limit  where  Im(k)  -*•  0.  For  the  purpose  of  numerical  integration, 
two  routes  present  themselves.  We  can  either  allow  k  and  a  to  have 
small  but  finite  imaginary  parts,  or  we  can  attempt  to  integrate 
in  a  sj^metric  fashion  about  the  singularities  at  z=a  and  z=k,  when 
they  lie  on  the  path  of  integration,  and  add  in  the  appropriate 
residues.  In  fact,  the  first  alternative  was  chosen.  It  was  evident 
that  the  factorization  functions  must  be  smooth,  hence  slowly  varying 
as  functions  of  complex  k  and  a.  Also,  experimentation  with  values 
of  loss  tangent  indicated  the  results  were  relatively  insensitive 
to  the  loss  tangent.  Hence,  for  purposes  of  numerical  integration 
complex  values  of  k  and  a  (being  k  and  o  )  are  defined  by 


kg  =  (1  +  i.005)k  (A-10) 

=  a  +  i.0025  k  (A-11) 

The  results  of  numerical  integration  for  values  of  n  (the  order  of 
the  Bessel  functions  listed  as  n  in  Table  A-1)  ranging  from  0  to 
3  are  presented  in  Table  A-1.  It  was  found,  by  way  of  confirmation, 
that  there  was  excellent  agreement  with  both  DC  and  asymptotic  forms. 
The  $M^  function  is  defined  by 


$M^(a)  =  a(k  +  a)M^(a) 


(A-12) 


It  remains  finite  as  k  0. 
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Table  A-1 

VALUES  OF  FACTORIZATION  FUNCTIONS  L+  AND  M+  COMPUTED  BY 
NUMERICAL  INTEGRATION 
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A-2.  Low-Frequency  Approximations 

Derivation  of  'DC  forms  for  the  factorization  functions  is 
not  difficult.  They  can  be  obtained  directly  from  the  small  argument 
approximations  for  Bessel  and  Hankel  functions.  For  n  ^  0,  we  get 
(Abramowitz,  Equations  (9.1.7)  and  (9.1.9)) 


L+(oi)L+(“a) 

=  iriJ^(Ya)Hj,^)(Ya) 

(A-13a) 

r^(f) 

(A-13b) 

.r  1 
n 

(A-13c) 

L+(a) 

✓n 

(A-14) 

M+(a)M^(-a) 

=  irij;(Ya)H^’)'(Ya) 

(A-15a) 

1  a 

i  r 

(A-15b) 

For  n  =  0,  we  note  that 

J^(2)  =  -  Ji(z) 

h[,^^'(z)  =  -  Hj^^(z) 


(A-lBc) 

(A-16) 


(A-17a) 

(A-17b) 
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Hence, 

M+(n=0,k,a)  =  L^(n=1,k,a)  (A-18) 

From  Abramowitz  Equations  (9.1.7)  and  (9.1.8),  for  n  =  0 

L+(ot)L^(-a)  =  fiiJQ(Ya)H^^^Ya)  (A-19a) 

^  TTi.l.  il  iln(Ya)  -  2jtn  ya  (A-19b) 

-ln(a^(k^-a^))  (A-19c) 

The  factorization  of  this  DC  form  is  non-trivial.  Based  onTWeinstein 

(1969),  p.  335-36]  we  may  approximate  L.(a)  ^  i/2  —  a/2k(kj^ 
a  «  k  «  1 .  /iln(2ka) 

Having  thus  obtained  the  first  low-frequency  approximation 
with  ease,  it  is  extraordinarily  difficult  to  obtain  a  better  approxi¬ 
mation  analytically.  Lee,  Jamnejad,  and  Mittra  present  a  simple 
derivation.  Unfortunately  it  does  not  agree  with  the  results  of 
numerical  integration.  This  can  be  seen  by  comparison  of  their 
formulas  with  the  plot  in  Figure  A-1,  which  shows  the  trajectories 
of  L^(k)  and  M^(k)  for  n=l  computed  by  numerical  integration.  The 
inclusion  of  additional  terms  from  the  small  argument  approximation 
does  not  improve  matters.  The  approach  taken,  therefore,  is  that 
for  small  values  of  ka,  the  values  of  and  are  computed  by 
numerical  integration,  and  interpolation  is  used. 

It  turns  out  that,  for  L^(n=l,o=k),  the  values  over  a  surprisingly 
large  range  (0  <  ka  <  2)  can  be  approximated  by 


L 


+ 


l-i.0584  ka 
(l+i.68  ka) 


(A-20) 


However,  this  does  not  admit  any  generalization,  nor  does  it  possess 
any  theoretical  justification. 


High  Frequencj 


3totic  Approximations 


Using  the  large  argument  approximations  for  Bessel  and  Hankel 
functions,  we  obtain  (from  Abramowitz  Equations  (9.2.17)  through 
(9.2.20)) 
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r 


F  igure  A-1.  Factorization  functions  L+,  M+,  $M+  for  n=1,a=k. 
Arrows  indicate  increasing  k. 
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(A-21a) 


L+(“)L+(-a)  =  ’fiJn(Ya)H^^)(Ya) 


i0n(Ya) 

=TriM^(Ya)cos8^(Ya)M^(Ya)e 

(A-21b) 

TTi  2  i20n(Ya) 

=  •^(M^(Ya))^(l+e  "  ) 

(A-21c) 

=  TTiJ’(Ya)H^^^'(Ya) 

(A-22a) 

i<l>n(Ya) 

=  TTiN^(Ya)cos4i^(Ya)N^(Ya)e 

(A-22b) 

TTi  2  i2<J>_(Ya) 

=  ■P-(N„(Ya))^(l+e  "  )  . 

(A-22c) 

Where,  from  Abramowitz,  Equations  (9.2.28)  through  (9.2.31) 


1  + 


4n2. 


8z 


0n(2)  <^2  -  +  + 


4n^-3 

8z2 


V^) 


(A-23) 

(A.24) 

(A-25) 

(A-26) 


For  L^(a)  we  observe  that  we  can  separately  factor  the  two  expres¬ 
sions,  making  L^(a)  the  product  of  two  factors,  both  of  which  are 
analytic  in  the  upper  half-plane. 


L+(a)  =  Li(ot)Lj^(a) 


(A-27) 
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(A-31) 


Mj(a) 


a(k+oi) 


■'^-  +  aa 


a(k+Qt) 


In  evaluating  and  M^^(a),  we  follow  the  derivation  of 

et  al,  (1975)  . 


Lj^(a)Lj^(-a)  =  1  +  e  " 


126  (Ya) 


M”(a)M|^(-a)  =  1  +  e  " 


i2<(>„(Ya) 


II- 


The  exact  expression  for  (a)  is  given  by 


Li^(o) 


=  exp< 


1  r  £n 


1  +  e 


i20n(Ya) 


z-a 


-  dz 


Let 


, ,  j 


z-a 


dz 


Then  the  identity 


ln(l  +  X)  =  ^  ;  lx|  <  1 

tn=l 


can  be  used  to  obtain 


i2me„(Ya) 

®  /  ixiiH-l  ®  „  n' '  ' 

1  =  1. ' 


m=l 


z-a 


dz 


provided  lm(0^)  >  0. 


Mittra 

(A-32a) 

(A-32b) 

(A-33) 

{A-34) 

(A-35) 

(A-36) 
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Let 


I 


m 


= ; 


i2m6  (ya) 
e  "  dz 


z-a 


(A-37) 


where 


e„(Ya) 

2 

-  /  n  .  1  \  ,  4n  •  1 

Ya  -  (^  +  7)11  +  — 

(A-38) 

ya  =  a' 

(A-39) 

z  =  k 

sinx 

(A-40a) 

dz  =  k 

cosx  dx 

(A-40b) 

ya  =  a^ 

/  2  2  2 

vk  -k  sin  x  =  ka  cosx 

(A-40c) 

The  contours  of  integration  in  the  z  and  x  plane  are  shown  in  Figures 
A-2  and  A-3,  respectively.  The  contour  in  the  x  plane  is  selected 
so  that  ya  has  the  right  sign,  and  the  integrand  vanishes  at  both 
ends  of  the  contour.  The  contour  in  the  x  plane  is  then  deformed 
to  the  steepest  descent  contour.  Since  no  poles  are  passed  over 
in  this  deformation  (for  Re(a)>0)  the  integral  is  unchanged.  Hence 


I 


m 


; 

c 


i2m0^{kacosT) 

ksinT-a 


kCOST  dT 


= ; 

^SOP 


cosx 

sinx-f 


i2me  (kacosx) 
e  "  dx 


With  9^(kacosT)  =  kacosx  -  (-j  +  + 


4n^-l 

Skacosx 


(A-41) 


(A-42) 


.  4n^- 
^  8(ka)^ 


(A-43) 
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Ijjj  is  an  integral  in  the  form 

f  F(t)  dT 

^SDP 


where  F(t)  =  — 

sinx 


a 

k 


K  =  ka 
f(T)  = 


i2m0^(kacosT) 

fa 


=  i2m(cosT  -  (7 


+ 

COST 


(A-44a) 

(A-44b) 


(A-44c) 


This  can  be  readily  evaluated  by  asymptotic  techniques.  We  make 
a  first  approximation  by  ignoring  the  pole  near  the  saddle  point. 


X,  sinx 

f'(x.)=0  =  i2m  (-sinx  + - 5 — ) 

^  cos  X 


(A-45) 


=  0 


(A-46) 


Then  f'Cxj)  =  i2m(-cosx+  ^'^3^  ) 


cos  X 


=  -i2m{l-XL) 


(A-47) 


The  asymptotic  value  of  the  integral  is  then 


<f(Ts) 


yicTTTj 


F(Ts) 


where  ^  is  the  angle  of  the  C^pp  through  the  saddle  point. 


(A-48) 


i2me„(ka) 


F(ts)  -  -  ^ 


f'Ctj)  =  -i2m(l-XL) 


(A-49a) 


(A-49b) 


(A-49c) 


Thence 


6 

m 


i2m0^(ka)g-'''^/'‘ 


(-  -) 
a 


(A-50) 


m-1  (Tl  'Ot' 

m=  I 


k  „ 
.r  —  e 
a 


-  i  7r/4 


Ji 


L-ir"  e 


i2in0„(ka) 


(A-51) 


The  indicated  sum  is  shown  by  [Chuang,  Liang,  and  Lee  (1975),  p.  773] 
to  be  a  modified  Lerch  function.  They  present  a  development  of  its 
properties  and  a  transformation  to  simplify  its  evaluation. 


r  (-ir  ~ 

— — ^75 -  =  L  - T77 


m=  I  m 


(A-52) 


where  +irmi  +  i2m0^(ka)  =  i2Tmip 


or  p  = 


6n(ka) 


(A-53) 


This  function  has  two  useful  and  simple  properties.  Since  the 
parameter  p  appears  only  in  the  exponent,  and 
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gi2iimp+i2iip  _  gi2iimp 
L(p+1.3/2)  =  L(p.3/2) 


(A-54) 


Furthermore 

L(l-p,3/2)  =  L*(p,3/2),  where  *  denotes  the  complex  conjugate. 

(A-55) 

By  using  these  relations,  it  is  possible  to  transform  any  value  of 
p  to  the  range  0  <  p  <  0.5.  In  this  range,  the  series  solution  is 
extremely  accurate 


L(p,v) 


^0  S'’ 


(A-56) 


For  \f3/2,  this  gives 


L(p,3/2)  =  -2Tr  »^.(l-i)  +  I  (R  +  il.p)p^'’  (A-57) 

n=0  "  " 


where  values  of  Rp  and  I  are  given  by  Table  I  of  Chuang,  Liang, 
and  Lee  (Note  sign  error^in  Equation  (35)  of  Chuang  et  al  (1975) 
corrected  here). 


We  can  similarly  find  an  asymptotic  form  for  M  (a).  The  analysis 

2 

is  identical  with  0„  (ka)  replaced  by  4)„(ka)  and  X. 

\(ka)  "  ”  8(ka)^ 

This  results  in  p  =  - - +  1/2.  Otherwise  the  expression  is 

unchanged.  ^ 


■  Li’(a)  ■ 

^L,M 

,  exp  ^ 

1  J 

(A-58) 


Where 


I 


L 


k  -iTi/4 
a  ® 


kin^ 


L(p,3/2) 


0n(ka) 

p  *  —  +  1/2  (A-59a) 
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■h  ■  5  /hoV  ^ 

A  comparison  of  these  results  with  those  obtained  by  numerical  inte¬ 
gration  is  instructive.  For  a=k,  the  agreement  is  very  good  down 
to  the  first  zero  of  J^(ka)  or  J^(ka)  (for  and  respectively) 

if  the  proper  form  of  ^^(ka)  or  <l>p(ka)  is  chosen.  This  is  illustrated 

in  Figure  A-4  for  M+(n=l,  a=k);  the  exact  value  of  is  compared 
to  that  obtained  by  using  the  two  term  and  three  term  approximations 
to  (|i,  (ka).  It  can  be  seen  that  a  rather  sharp  minimum  occurs  for 
the  4xact  form  of  M_|^  for  that  value  of  ka  for  which  J,'(ka)=0.  For 
the  two-term  and  three-term  forms  of  (|)i(ka),  the  minimum  occurs  at 
the  value  for  which  cos({»i (ka)=0.  Thus,  the  more  accurate  9p  or  <|)p 
should,  in  general,  lead  to  a  more  accurate  approximation  for  L. 
or  M^. 


and  are  also  functions  of  ot/k;  we  must  explore  the  ac 
curacy  of  this  approximation  when  this  ratio  varies.  This  is 
tant  since  a/<  =  cos9  for  scattering,  radiation,  and  coupling 
.  ^ 

and  (Vk 


impor- 

problems. 


)  for  coupling  and  radiation  of  waveguide  modes. 


In  Figure  A-£,  the  value  of  L^(n=l)  is  plotted  in  the  complex  plane, 
where  the  main  curve  traces  out  the  real  and  imaginary  parts  of 
with  o^k  and  Q/\  as  a  parameter.  The  dashed  curves  then  indicate 
the  tr^ectories  which  occur  when  varying  a/k  while  holding  ka  (or 
D/X)  fixed.  It  is  immediately  evident  that  the  preceding  approxi¬ 
mation  does  not,  in  general,  capture  the  nature  of  this  behavior. 

The  reason  for  this  can  be  seen  from  Figure  A-3.  As  a/k  varies, 
the  pole  located  at  Xp  =  sin  (a/k)  moves  closer  to  the  saddle  point 
at  15=0.  We  totally  ignored  this  pole  originally.  Without  belaboring 
the  point,  it  will  only  be  stated  that  rederiving  the  approximation 
incorporating  the  effect  of  the  pole  produced  only  a  marginal  im¬ 
provement.  It  was  therefore  necessary  to  include  two  terms  in  the 
asymptotic  expansion.  This  derivation  follows. 


With  the  same  definitions  of  <,f(T)  and  F(t),  we  have  the 
asymptotic  approximation 
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Figure  A-4.  Factorization  function  M+  for  n=l,  a=k 
calculated  by  numerical  integration  and  asymptotic 
approximation  based  on  two  and  three  term 
expressions  for  phase  of 
Bessel  function. 


Factorization  function  L+  for  n=l.  Solid  curve  gives  a=l<,  varying  k 
curves  vary  a/k  holding  D/X  fixed  to  the  values  specified  by  arrows. 

□  =  exact  value  (by  numerical  integration) 

X  =  first  asymptotic  approximation. 


icf(T  )  i(l) 

I  e  ^  e 
m 


2iT _ 

<7^ 


F(t3) 


+  iA  (1  -  f(<A)) 


I 


+  F'(Tg) 


f'"(Tj 


{A-60) 


where  f(X)  is  the  transition  function  defined  and  discussed  by 

[Kouyoumjian  and  Pathak  (1974),  p.  1453].  A  is  a  distance  parameter, 

relating  to  the  separation  of  x  and  x  ,  defined  by 

s  p 


A  =  i(f(x^)  -  f(Xp))  (A-61) 

The  saddle  point  is  unchanged,  so  Xg=0.  The  pole  is  located  at 

Tp  .  s1n-'  2 

We  thus  derive  Table  A-2. 


(A.62) 


TABLE  A-2 

SUMMARY  OF  ASYMPTOTIC  FUNCTIONS  AT  SADDLE  POINT  IN  t-PLANE 


Function  form 

Value  at 

f(T) 

i2m(cosT-(^  + 

i2m(l-(^  + 

f'(T) 

i2m(-sinT+X,  ) 

^  COS'^T 

0 

f(T) 

2 

i2m(-cosT+X.  -1-  -^) 

COS'^T 

i2m(-1+X|^) 

f"  (t) 

i2m(sinT+X.  -^^”1  (6-cos^t)) 

^  COS^T 

0 

12m(cosT*x,  ':<'5S-20co52t*24) 

^  COS^T 

i2m(l+5XL) 

F(t) 

COST 

sinx-a/k 

-k/a 

F'  (t) 

a/k  sinx-l 
(sinx-a/k)^ 

ro 

F"{t) 

2 

-a/k  sinxcosx-(a/k)  COSX+2COSX 

^  -  2*^ 
a  'a 

(sinx-cx/k)^ 

-V  k  /  ir 

a/ITan^ 


,-in/4  „-I/2 


1 


1  1+5Xl  1-2(i) 


2. 


L**  (-HXl)‘  -HXl  J 


=  wS 


=  2 


/ 


l-(a/k)^  -  1  +  \(  ^ 


.yr 


(o/k)' 


-3+9Xl  +  8(^)  (l-X^) 

8(1-Xl)^ 


_  k  /  if 
-  a/TrarT^ 


,-iTr/4 


A,  i2me_(ka)  f  ] 

>r - -e  +  AQ6(l-lf(mka6))> 

i/fiT  k  J 


00  f  tvITH-I 

=  I 


m=l 


m  m 


^  m= 


»r  A 


00  .  .im-l  (-A,)  i2nie  (ka)  f  'j 


®  (  ,,m  i2me  (ka)  f  'I 

1  ^  -^TT®  <l+Ap6(l-f(mka6)) 

m=  1  \  J 

^  A^  |L(p,3/2)+Ag6G(p,q)j 


(A-65b) 

(A-66) 

(A-67) 

{A-68) 

(A-69) 


(A-70) 
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where 


G(p.q)  =  ?  m-3/2v2™P(l-f(mq))  (A-71) 

m=1 

0n(ka) 

p  =  -  +  1/2  (A-72a) 

q  =  ka(S  (A-72b) 

We  could  follow  exactly  the  same  analysis  to  obtain  M+,  except  that 
0n(l<a)  -*■  4)^(ka) 


m 


4n^+3 

8(ka)^ 


(A-73) 


These  formulas  produce  a  very  good  agreement  with  the  results  of 
numerical  integration,  down  to  cx/k=0.2. 


s 


-  2i 

4)(|)  ”  k(1+cose^)(l+cos0^) 


Jq  epcos 


J^(kasin0.) 

H_j^(kcos0^.) 


J^(kasin05) 

M^(kcoseg) 


(l+cos0^.)(l+cos0g) 

2(cos0.+cos0g)  i_^2 

'  n  - 


(B-4) 


^0({)  k(l+cos0^)  "‘'^s  s1n0gL^{kcos0g)  M^(kcos0^)  * 

n 

(B-5) 


.  41  ? 

^4>0  k( l+cos0j)  M^(kcos0g)  sin0^.L^(kcos0^-) 

n 

{B-6) 

Where 

<|)^  =  TT  is  assumed 
nL^(k) 

^  “  ?kiW^(k) 


1 


e 


n 


n  =  0 

n  =  1,2,3... 


(B-8) 


It  was  found  that  Chuang  et  al  (1975)  contained  several  misprints. 
A  correction  letter  is  reproduced  on  the  following  page. 
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AP1M4 


Cal.  1 


Coirecdon  to  “High  Frequency  Scattering  from  an  Open- 
Ended  Semi-Infinite  Cylinder" 

C.  A.  CHUANG.  MEM  BEK,  IEEE,  CHARLES  S.  LIANG, 
MEMBER,  IEEE,  AND  SHUNG-WU  LEE, 

SENIOR  MEMBER.  IEEE 

In  the  above  paper* ,  there  were  six  sign  misprints. 

1)  Time  factor  should  have  read  exp(  -/tor). 

2)  In  (1 ),  X  should  have  read  (-jr). 

3)  In  (4),  (1  +  cos  d)  in  numerator  should  have  read  (1  - 
cos  0). 

4)  In  (7),  (cos  0  +  cos  0o)  should  have  read  (-l)(cos  0  + 
cos  0o).  ■ 

5)  In  (1 1),  (a  -  Of')  should  have  read  (a  -  a). 

6)  In  (35),  (I  +  i)  should  have  read  (I  -  i). 

the  above  sign  misprints  do  not  affect  any  other 
equations  or  numerical  results. 

The  authors  wish  to  thank  T.  W.  Johnson  and  D.  L.  Moffatt 
for  bringing  some  of  these  errors  to  their  attention. 

C.  A.  Chuang  is  with  Aeionutronic-Ford/i^’DL,  Palo  Alto,  CA.  He 
is  now  at  552  Solitaire  Dr.,  Indialantic,  FL  32903. 

C.  S.  Liang  is  w  ith  the  Fort  Worth  Division,  General  Dynamics,  Fort 
Worth,  TX  76101. 

S.  W.  Lee  is  with  the  Department  of  Electrical  Engineering,  Univer¬ 
sity  of  Illinois,  I’rbana,  IL  61801. 

•  C.  A.  Chuang,  C.  S.  Liang,  and  S.  W.  Lee,  /fff  Trans.  Antenna 
Propagat.,  vol.  AP-23,  pp.  770-776,  Nov.  1975. 
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B-2.  Coupling  of  Incident  Field  to  Waveguide  Modes 


With  the  same  definition  of  the  incident  field  (as  in  Equation 
(B-2))  the  axial  fields  within  the  waveguide  can  be  represented  as 
follows; 


E 


z 


00 


I 

m=l 


A^„sin  n4>  i  '^rr^niTr 

I  nm  j 


(B-9) 


00  00 
1 

n=0  m= 


*  A 

B„„sin  n4> 
nm 

B^  cos  n()) 
nm 


•■'“nmZ 

nm 


(B-10) 


'^n^-^nm 


)=0 


n'‘'nm' _ 


'nm 


nm 


The  top  line  of  Equation  (B-9)  indicates  that  Eq  produces 
E^  with  cos  n(|)  dependence 

The  bottom  line  of  Equation  (B-9)  indicates  that  e!  produces 
E^  with  sin  r\^  dependence  ^ 


e  2e_i'^j„  L.(a  )  J  (kasinS.) 
I  -  n  nm  nm  n  _  .  ^ 

nm  ~  "  a„„a  L (kcose .)  kasine . 

nm  +  1  1 


77?  Z(o.„„-kcos6,) 
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(B-I2) 


A*  .  -  J™ 

nm 


-j;!" -  j;(k»sine,)  ^ 

ka  a^(l+cos6.)  M+^^cose.)  n  ^ 

nm  '  V  r 


„e  , 

“nm  °  * - 'J~ 

■^nm 


M4.(“J,m}  J_(kasin0.)  f„ 
+  _  nm  n^  i '  n 

L^(kcos6^)  kasinS^. 


(B-13) 


J***  = 
nm 


2e„i^(k+a'  ) 
n  '  nm' 


^)(Hcose.) 

■^nm 


”+(“nm^ 

H^(k^s6,)  J;(k«ine,)  . 


fn  (lc.a;^)(Hcose.) 
Nf!  2(kcose,-a;^) 


(B-14) 


B“3.  Radiation  Patterns  From  an  Open  Waveguide 

The  incident  field  is  assumed  to  be  a  single  waveguide  mode, 
either  TE  or  TM,  with  axial  fields  given  by 


.i  _  ,  _ '“nm' 

n'  nm' 


(B-15) 


Jn(jLP/3)  iot*  z 

Z  '  J7j'  ) - ® 

n  nm' 


(B-16) 


The  radiated  field  is  given  by 
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(B-17) 


^  _  (-D^^^e^k^^cos  n*}.^  J„(kasin95) 

^0E  t^(kcos9g)  kasinSg 


2(%,-l<cosSj) 


(B-18) 


eH 


(-^"■^^^a^sin  n(j)g  J^(kasineg) 

(k-a^)  L^Ckcose^)  kasinO^  ^[7^ 


(B-19) 


(-i)'’*'a^ksin  L^(cy„) 


j^(Ucose,)  J;(kas1ne^)  ^  (B-BO) 


'<PH 


(-t)"*'e„kacos  0*5  M^(ay  . ^ 

2(i+cos05)(a;i^-k)  F^Tkcoie^)  * 


i  .  (K-c^X'^cose^) 

77?  Z(kCOS03-OC^„) 


(B-21) 
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B-4.  Internal  Reflection  and  Mode  Conversion 


field 


With  incident  field  as  described  in  Section  B-3,  the  scattered 
is  given  by 


s  _  . 


.-'V 


"III”"  "♦  C"’" 


"nm 


nm 


(B-22) 


Is 


„s  .  r 

■  I  " 


Tjj^sin  n4)  rJJJcos 


'mt" 


nm 

H 

11  S 


(B-23) 


'  ■  S;: 


ynH  _  ^ng,‘^**^+^°nm^‘'+^°Vn!,^  ^n 

1-ff 


l-f?  2k(a^„+a^^) 


(B-24) 

(B-25) 


„E  ka(k+a;^) 


m£ 


2  \  *-+^®nm^*^+^“nJl^  ,  -2 


‘n£ 


I-  2 
'nl 


1-f! 


(B-26) 


mi 


°‘nl^''"“r»nM^  ~  2 

^nl 


(B-27) 
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APPENDIX  C 

ELEMENTS  OF  DYADIC  GREEN’S  FUNCTION 
FOR  A  CIRCULAR  WAVEGUIDE 


The  dyadic  electric  Green's  function,  omitting  the  source- 
singular  term,  is  given  by  Tai  (1971).  It  is  assumed,  in  this  appen¬ 
dix,  that  the  azimuthal  dependence  for  Ep,  E  ,  Op,  J  ,  and  H*  is 
cos  n<|),  and  for  Ea,  J*,  Hp,  and  H  is  sin  n(|).  This  azimuthal  depen¬ 
dence  is  understood  and  suppressed.  The  summation  over  n  is  also 
suppressed. 
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